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Abstract 



We describe an analytic continuation of the Euclidean Grosse-Wulkenhaar and LSZ models 
which defines a one-parameter family of duality covariant noncommutative field theories inter- 
, polating between Euclidean and Minkowski space versions of these models, and provides an 

' alternative regularization to the usual Feynman prescription. This regularization allows for a 

\^ , matrix model representation of the field theories in terms of a complex generalization of the 

usual basis of Landau wavef unctions. The corresponding propagators are calculated and iden- 
, tified with the Feynman propagators of the field theories. The regulated quantum field theories 

\C • are shown to be UV/IR-duality covariant. We study the asymptotics of the regularized propa- 

, gators in position and matrix space representations, and confirm that they gcnerically possess 

a comparably good decay behaviour as in the Euclidean case. 



^ ■ 1 Introduction and Summary 



This paper is devoted to an in-depth study of the perturbative properties and renormalizability of 
noncommutative (p*'^-type scalar field theories on real vector spaces subjected to a Moyal deforma- 
tion. The vast majority of the hterature on this subject has been devoted to Euclidean quantum 
field theory. The most prominent feature of these models is the notorious mixing of ultraviolet and 
infrared modes, which renders the noncommutative (j)*^ field theories nonrenormalizable [M"VRSOO] 
(see e.g. [Sza03] for a review). Grosse and Wulkenhaar demonstrated how to obtain field theories 
which are renormalizable to all orders in perturbation theory by extending the kinetic term of 
the (j)*'^ Lagrangian by an additional harmonic oscillator potential [GW03, GW05b]. The Grosse- 
'Wulkenhaar model also has vanishing beta-functions and its perturbation series is likely to be 
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Borel summable [Riv07]; in two dimensions this has been recently estabhshed in [Wanll]. In four 
dimensions, the Euchdean Grosse-Wulkenhaar model is the first rigorous four-dimensional quan- 
tum field theory without unnatural cutoff which is expected to exist non-perturbatively and is not 
asymptotically free. 

The continuation of the Euclidean models to noncommutative Minkowski space is presently an 
open problem which is plagued by both conceptual and technical difficulties. The original problems 
were unveiled in [GMOO], where it was found that the standard perturbative expansion in terms 
of Feynman diagrams leads to a violation of unitarity if space and time do not commute. As 
subsequently pointed out in [BDFP02], this due to the failure of Wick's theorem, which does not 
apply to non-local interactions in general. By using canonical quantization in the Hamiltonian 
framework involving the Dyson series and time-ordered products of the interaction Hamiltonian, 
the resulting field theory is still unitary but no longer equivalent to the Lagrangian formulation 
of the quantum field theory in the path integral framework. For models built on the Hamiltonian 
framework see e.g. [DFR95, Bah04, PialO]. Yet another inequivalent perturbative approach is based 
on the Yang-Feldman formalism [BDFP02, Bah04], which also gives a unitary noncommutative 
quantum field theory on Minkowski space with time-like noncommutativity.^ 

The UV/IR-mixing problem of the ordinary (/>*"" field theory is absent in the Hamiltonian 
framework to lowest orders, and it has long been an open question as to whether it exists at all. 
Only recently has UV/IR-mixing been shown to still occur, albeit through a mechanism which is 
different from that of the Euclidean setting with modified Feynman rules [BahlOb]. It has also 
been shown that UV/IR-mixing arises in the Yang-Feldman formalism [Zahllb]. 

Since the perturbative setups in the Hamiltonian and Yang-Feldman formalisms are quite com- 
plicated, it would be desirable to have an equivalent Euclidean path integral formalism which 
simplifies the combinatorial aspects of perturbation theory. However, the relationship between 
the Euclidean and Minkowski space theories when time and space do not commute is unclear. 
In [BahlOa] it has been shown that the Euclidean counterparts of the n-point functions for the 
Klein-Gordon theory on noncommutative Minkowski space in the Hamiltonian formalism are not 
those which follow from the standard Euclidean framework, but appear with on-shell twisting 
factors involving only on-shell momenta. 

In this paper we will pursue the other direction of this correspondence, starting with specific 
field theories in Euclidean space. In order to find noncommutative field theories in hyperbolic 
signature which are free from UV/IR-mixing, we construct Minkowski space counterparts of the 
Grosse-Wulkenhaar model and its generalizations known as the LSZ models [LSZ03, LSZ04]. These 
models differ from the standard noncommutative (p*'^ field theory by the introduction of an external 
background field into the kinetic term of the Lagrangian, making them covariant under the duality 
comprising Fourier transformation plus a rescaling of the fields [LS02]. This duality is believed 
to be related to the improved asymptotic behaviours of the propagators, which suppresses the 
UV/IR-mixing. This UV/IR duality may thus be responsible for the renormalizability of these 
field theories. The Euclidean Grosse-Wulkenhaar and LSZ models are defined via path integral 
quantization, which leads to a violation of unitarity for the usual noncommutative field theories in 
Minkowski space. Here we will be interested in the renormalization properties of their hyperbolic 
counterparts; unitarity of these quantum field theories will be addressed elsewhere. 

In Euclidean space the introduction of an external field has the useful additional effect that 
the corresponding wave operators have discrete spectra and the models can be analysed with the 

^ Still another approach is provided by the twist-deformation formalism for noncommutative quantum field theory 
on Minkowski space which is considered in [BGM+07, GL07, GL08, S0IO8]; here one first quantizes the classical field 
theory before deforming spacetime, and the free part of the quantum field theory also differs from its commutative 
counterpart. 
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help of a matrix basis for expansion of fields; the countably infinite set of eigenfunctions are the 
Landau wavefunctions which diagonalize the free parts of the action. This basis defines a mapping 
of the duality covariant field theories onto matrix models, which permits a simple and natural 
regularization of the field theories while maintaining duality manifestly at quantum level. In this 
way Grosse and Wulkenhaar were able to prove the renormalizability of their model to all orders of 
perturbation theory. In addition, it has been used to solve the LSZ model exactly and demonstrate 
the vanishing of the beta-function. 

However, in passing to hyperbolic signature, the background field, which is a magnetic field in 
the Euclidean metric, now plays the role of an electric field. This yields a qualitative change due to 
the work done on the particles by the field. The electric field accelerates and splits virtual dipole 
pairs leading to pair production. This is reflected in the spectra of the wave operators, which now 
have a continuous part and are unbounded from below. 

In [Zahlla] the perturbative expansion in terms of modified Feynman diagrams in the contin- 
uous eigenvalue representation has been investigated. At one- loop order unusual divergences arise 
which are very likely to be non-renormalizable. Moreover, the retarded propagator in position 
space is no longer a tempered distribution in general. 

In [FS09] a different approach has been investigated, where a set of resonance states has been 
used to expand the field theory in a discrete set of functions. In the following we will take yet 
another path, which is related the resonance expansion found in [FS09], but which potentially 
avoids the associated technical problems. We will show that the Grosse- Wulkenhaar and LSZ 
models allow for well-defined analytic continuations to Minkowski space with the help of a special 
regularization, that we call the "■iJ-regularization" , which is a suitable replacement for Feynman's 
i e-prescription. 

This approach may also avoid the strange divergences found in [Zahlla]. These divergences 
come from squares of Dirac delta- functions which arise from undetermined loop integrations. They 
are not ultraviolet divergences in the usual sense, as they occur before performing loop integrals, 
and they show up in every (/)*"-theory with n > 3 for graphs with an unbroken internal line. Using 
the 'i?-regularization instead of the usual ie regularization, one gets a discrete spectrum instead 
of a continuous spectrum, leading to Kronecker delta-functions and sums rather than Dirac delta- 
distributions and integrals. This procedure renders these diagrams finite, and at the same time 
keeps the model duality covariant. 

The outline of this paper is as follows. We will define duality covariant quantum field theories 
on Minkowski space based on the work [FS09, FSll], and describe some of their renormalization 
properties. In order to employ an expansion of the action functionals in terms of the resonance 
states found in [FS09], we will regularize the models such that the resonances turn into genuine 
eigenfunctions of the regularized wave operators; this new matrix basis and the corresponding 
matrix model representations of the LSZ and Grosse- Wulkenhaar models on Minkowski space are 
described in detail in §2-§3. These wave operators are related by Weyl-Wigner correspondence to 
the complex harmonic oscillator Hamiltonian, which interpolates between the ordinary and inverted 
harmonic oscillator Hamiltonians, and thus between the Euclidean and Minkowski space theories; 
this unifies both theories into a one-parameter family of duality covariant noncommutative quantum 
field theories. The Feynman graphs are analytic continuations of the Euclidean diagrams. We show 
that this regularized matrix basis is a bi-orthogonal system whose linear span is the space of square- 
integrable functions. At the quantum level and in the limit of vanishing electric background, this 
regularization turns into the usual i e-prescription. For the special case of a Klein-Gordon theory 
in a constant external electric field, where the various propagators are known, we recalculate in §4 
the propagator using the complex matrix basis and verify that the regularization leads to Feynman 
propagators. This confirms the equivalence to the i e-prescription, and demonstrates that the ■!?- 
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regularization is also connected to causality of the quantum field theory. Using the i?-regularization, 
we show in §5 that a cutoff can be introduced that renders the duality covariant field theories finite 
at every order of perturbation theory and at the same time imposes duality covariance manifestly. 
In §6 we derive the propagators for the regularized models which include the Euclidean space 
propagators and the Minkowski space causal propagators as special cases; away from the hyperbolic 
point our propagators have a good decay behaviour in all directions and are singular at coincident 
positions. The ??-regularization turns out to improve their asymptotic behaviour and may thus 
be crucial for the renormalization programme. However, due to the oscillatory behaviours of the 
occuring integrands in Minkowski space, the corresponding asymptotics are much more difficult to 
derive than in the Euclidean case. For a special case of the LSZ model we find that the exponential 
decay in the short Euclidean space variables ceases if one goes over to Minkowski space, but 
persists in a neighbourhood of the hyperbolic point in the one-parameter family of field theories. 
The regularization thus gives a means to control the decay behaviour of the propagators. The 
applicability of the matrix basis in this context, however, is still an open question; the detailed 
analysis of the removal of the matrix regularization will not be addressed in this paper. As we discuss 
in the following, the Minkowski limit of our one-parameter family is very singular; see [Fisll] for 
a detailed analysis of some of the uncontrollable divergences which arise. In the following we will 
simply regard the i?-regularized field theories as the appropriate well-defined analytic continuations 
of the Euclidean space models. 

The derivations of propagators with the help of the matrix basis may be compared to calculations 
using other methods, such as Schwinger's proper time formalism [SchSl], the "sum over solutions 
method" [FGS91], or the eigenvalue method using the continuous eigenbasis [Rit78]. Compared to 
the latter technique the matrix basis involves only polynomials and sums instead of complicated 
integral expressions, and thus brings along a huge simplification. In §4.2 the causal propagator for a 
massive complex scalar field in four dimensions in the background of a constant electric field in the 
"iJ-regularization is computed. As a further demonstration of how the matrix basis can be applied, we 
calculate the one-loop effective action of the Klein-Gordon theory in a background electric field (see 
Appendix B). Finally, in Lemma D.13 we demonstrate that the -i^-regularization can likewise be used 
regulate the standard mass-shell singularities in the Feynman propagator for the free Klein-Gordon 
theory. We propose that going beyond the case of a constant background field might be possible 
using our alternative regularization and the matrix basis, by perturbing varying field configurations 
around a uniform background. This might help in probing quantum electrodynamics in the non- 
perturbative regime (see e.g. [Rin03, III09, Dun09, ILMIO]). We conclude that the matrix basis 
may serve as a powerful computational tool in simplifying some otherwise cumbersome calculations. 

2 Covariant Relativistic Noncommutative Field Theory 

In this section we will introduce the duality covariant models in Minkowski space. We show that it 
is possible to construct a well-defined matrix model representation of the corresponding quantum 
field theories through a suitable regularization, that we call -regularization, which is an alternative 
to the usual i e-prescription. For this, we will use the Weyl-Wigner transformation to map the 
eigenvalue problem of the ^^-regularized wave operators to that of the complex harmonic oscillator. 

2.1 Formulation of the Duality Covariant Models 

We work in D = 2n spacetime dimensions with metric of signature (1,— 1,...,— 1). Tensors will 
be labelled by Greek indices /x, i^, . . . ranging from to d = D — 1. Throughout we use the Einstein 
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summation convention. For simplicity we will denote the hyperbolic norm square of vectors a = (a'^) 
as 

II Il2 2 2 2 a 2 /n i \ 

l|a|lM = «o-ai = a^a^" =: . (2.1) 

Euclidean space dimensions are labelled by Latin indices i, j, . . . ranging from 1 to D, and norm 
squares of vectors a = (a*) with respect to the D-dimensional Euclidean metric are denoted 

l|a|lE = a?H bal) = Qia' =: . (2.2) 

Position vectors are denoted x = (x^), with derivatives := ^|jr; in two dimensions we often 



write X = {t,x). The dual pairing between a covariant vector x = (x^) G M and a contravariant 
vector k = (A;^) G (IR-^)* is written k ■ x = x'^. 

The LSZ model is a complex </)^-theory defined by the action 

'5lsz = So + 5int (2.3) 

with 



So 



d^xc^*{x)(aKl + {l-a)k^ 



(^) , 



d^x ( a {(p* ★e (j) -ke cf)* *e </>) (^e) + /3 (</>* *0 (/>* ★e *e </>) {x) 



(2.4) 



(2.5) 



where o" G [0, 1], a, /3 G M+ := [0, oo), and > are the mass and coupling parameters. The 

generalized momentum operators Km and generalized dual momentum operators are given by 



i9„ 



and they obey the commutation relations 



and 



K, 



-2iF, 



and 



[K^, K,] = . 



(2.6) 
(2.7) 



The star-product of arbitrary Schwartz functions f{x),g{x) G 5(M^) is given by 



(/*e g)ix) :-- 



1 



vr^ldetei 



j d'^y I d''zf{x + y)g{x + z)e-^'y'^''^ = {g**en{x) (2.8) 

with respect to a constant, real-valued, antisymmetric and non-degenerate D x D deformation 
matrix Q. 

The coordinate system is chosen such that G takes the canonical skew-diagonal form 

/ 00 \ 



G = (G 



-Oo 



Oi 
-Oi 



V 



On^l 



(2.9) 



with 6k > for k = 0, 1, . . . , — 1. The constant electromagnetic field strength tensor is 
likewise given by 



/ E 

-E 



\ 



Bi 
-Bi 



V 





-Bn-i / 



(2.10) 
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with E, Bk > and 

Eeo = Bkek = 2n (2.11) 

for = 1, . . . , y — 1 and < 17 < 1. We will sometimes regard and F as invertible linear 
maps e : (M^)* R° and F : ^ (M^)*. The LSZ models with a = I, where only the 
generalized momentum operator appears, are called "critical" models, while those with = 1, 
where the field theory is invariant under the UV/IR duality between position and momentum space 
representations, are called "self-dual" models. 

We will solve the eigenvalue equation for the wave operator and relate it to that of the Euclidean 
case. For this, we note that the 2n-dimensional wave operators break up into n blocks with 

n— 1 n— 1 

K' = E(P')fc and K2 = j;(p2), . (2.12) 

A:=0 fc=0 

The operators 

(P')fc = (dl + dik+i) + 2iBk {x''+' d2k - d2k+i) - Bl (xlk + xL^J , 

(P')fc = + ^Ifc+i) -1\Bk ix''+' - d,k+i) - Bl {xlk + xL+i) 

for k = 1, . . . ,n — 1 act on two-dimensional Euclidean Klein-Gordon fields in a constant external 
magnetic background of field strengths ±22?^, respectively, while the operators 

(Pu)o = -{dl - dl) -2iE{x'do+x' di) - E' {xl - xl) , 



(2.13) 



(P^')o = -{dl - Of) + 2iEix'do+x° dl) - E^ (xl 



(2.14) 



act on 1+1-dimensional Klein-Gordon fields in a constant electric background with field strengths 
± 2E, respectively. Since all component operators (P^)fc and (P^)a: for /c = 0, 1, . . . , n — 1 mutually 
commute, the diagonalization of the full wave operators amounts to diagonalizing each of their 
two-dimensional blocks independently. 

As is well-known the spectra of the operators (2.13) are discrete with corresponding eigenfunc- 
tions the Landau wavefunctions fmn satisfying 

{P^^)k f^n'Hxk) = -4Bk{m + ^)f^^^\xk) , 

{^l)kf!^^\xk) = -4Bk{n + \)f'^^\xk) (2.15) 

for m,n & No, where we write = [x"^^ , x"^^^^) G for /c = 0, 1, . . . , n — 1. These functions are 
Wigner transformations of tensor products of harmonic oscillator number basis states 

f^J{x,y)=\N[\m){n\]{x,y) (2.16) 

of frequency Q,, where, in two dimensions with deformation parameter = 9, the Wigner dis- 
tribution function of a compact operator p on Fock space is the Schwartz function on R? given 
by [Sza03] 

\N[p]{x,y) = j dk e'^y/^ {x + k/2\p\x-k/2) . (2.17) 

On the other hand, the spectra of the operators (2.14) are continuous, with corresponding eigen- 
functions given by Wigner transformations of tensor products of parabolic cylinder functions [FS09, 
Zahlla], denoted Xpq with g S R, which solve the eigenvalue equations 

{^l)oXpq{xo) = "^Epxpqixo) and {PI)q Xpqixo) = 4.E q Xpqixo) . (2.18) 



6 



For generic a the free part of the LSZ action (2.4) can be rewritten as 

5o = y d^x</)*(^) {Y.X^^^V?~xl-iJ?) <^{x) (2.19) 

with F = {2a — 1) F = {2a — 1) {F^^,,^) and x^j = 20~J^ x'^. The free action thus describes a massive 
complex scalar field coupled to a constant electromagnetic background and in an oscillator potential 
proportional to Q,"^ xj^. The Grosse-Wulkenhaar model in D = 2n spacetime dimensions is the LSZ 
model for £7 = ^ and a = /? = ^ with real scalar fields. The action is thus 

5gw = y d^x ^ <i){x) {-dl + n^S:l- /i') H^) - 9 J d^'x {4> *e *e *e 4>){x) ■ (2.20) 

The dimensional wave operator again reduces to a sum of n — 1 Euclidean wave operators plus 
a two-dimensional wave operator in Minkowski signature 

h (P')o + i (Pm')o - = -id'o - dl) - (xg - xl) - ^2 (2.21) 

with frequency = E9q/2. The main difference, besides the hyperbolic signature, is an extra 
minus sign in front of the O-term. The corresponding wave operator is given by the Hamiltonian 
of a harmonic oscillator with imaginary frequency, known as the inverted harmonic oscillator. 

The corresponding models on Euclidean space are defined in terms of the wave operator cr K? + 
(1 — o") K? + which also split up into n blocks made up of the operators 

n n 

K? = E(P?)'^- K? = E(P?)'^- (2.22) 

k=l k=l 

After relabelling of coordinates, one can relate {P'f)k = — (P^)fc-i and (Pf)fc = — (P^)fc-i for k = 
2, . . . ,n, whereas {Pf)i and {Pl)i are of the same form as (2.13). Thus in contrast to the mixed 
discrete and continuous spectrum of the hyperbolic space wave operator, the Euclidean case deals 
with purely discrete spectrum. This situation is responsible for the powerful application of the 
matrix model representation of Grosse and Wulkenhaar [GW03, GW05a, GW05b]. 

The duality covariant field theories involve two parameters and F. In the commutative limit 
= 0, one recovers the field theory for an interacting scalar field in a constant electromagnetic 
background; in §4.2 and Appendix B we demonstrate how to reproduce the known standard results 
in the literature using the novel regularization we propose below. In the vanishing background limit 
-F = 0, we recover the usual (/)^-theories on noncommutative Minkowski space together with their 
UV/IR mixing problems as discussed in §1; in Lemma D.13 we illustrate how our regularization 
is applicable in this case as well. Neither of these two limits possess duality covariance. In the 
self-dual limit F = (20)~^ the field theory is duality invariant; the matrix representation we obtain 
below at the self-dual point makes no sense in the limit F = Q. 

2.2 Spectral Decomposition and i?-Regularization 

The external electromagnetic background will be treated by considering all terms quadratic in 
the fields as being part of the free action. Then the path integral quantization gives the usual 
(modified) Feynman diagrams but with the dressed propagator for the scalar field moving in this 
background. It is a feature of most field theories defined on hyperbolic space that there is more 
than one propagator, i.e. a distribution whose kernel A(a;, y) solves the partial differential equation 
Da;A(ic, y) = 5{x — y) with Dx the wave operator of the field theory. This is due to the occurence of 
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zero eigenvalues of Dx, which prevents the naive inversion of the operator to give a propagator. It is 
therefore necessary to impose further conditions so as to make the solution of this problem unique. 
This may be done either by imposing boundary conditions, by postulating a spectral representation, 
or by extending the wave operator so as to make the solution of the partial differential equation 
unique. 

For the ordinary scalar field theory, the i e-prescription is a method to single out a specific 
propagator, namely the Feynman propagator. In the commutative field theories, this prescription 
enhances the action by an additional term \e J (jp' which for e > ensures the required asymp- 
totic damping of the integrand in the partition function at — t- oo (rather than an oscillatory 
behaviour) , and at the same time regularizes the singularity of the free propagator and furthermore 
imposes causality. In this particular case it is also the infinitesimal version of the Wick rotation to 
Euclidean space 1 1— )• e^'^t. 

However, in our case the field theories defined on the two different spacetimes are not related 
by this ordinary Wick rotation - it has to be accompanied by an additional transformation E i— >• 
zbiS. This is not surprising, since the model can be viewed as a field theory on a curved non- 
stationary spacetime, for which this is a generic feature [DeW75]. Another characteristic of those 
field theories is that the multitude of different equivalent definitions of the Feynman propagator is 
resolved [CR77]; we return to this point in §4.1. Since we are interested in an analytic continuation 
of the Euclidean space models in a path integral framework, it is the propagator we obtain by this 
transformation that we are concerned with. The extra transformation of the magnetic field strength 
is also in harmony with the fact that in order to ensure the commutation relation [x^ , a;*] = i 0*^* for 
both Euclidean and Minkowski space, the deformation parameter Q^'^ has to transform accordingly 
to compensate the phase coming from the Wick rotation. For the duality invariant field theories, 
i.e. at the self-dual point = 1, the deformation matrix is proportional to the field strength tensor, 
which in turn implies a rotation of the field strength. 

In [FS09] it was shown that the actions of the Minkowski and Euclidean space wave operators, 
(P^)o and (Pf)i, can be represented as a star-product with a classical Hamiltonian at the self-dual 
point = 1. To compare to the Euclidean version, we have to identify B = E and the ordered 
coordinate pairs x = {x^,x^) = {t,x) to find^ 



(P,?)i fix) = E^ [x^ + fix) and (P^2)o fix) = E^ (x^ - t^) ^^/^ /(^) (2.23) 



(P?)i fix) = fix) *2/E + t^) and (p2)o fix) = fix) " , (2-24) 



Consequently, there is a one-parameter family of operators which continuously interpolates between 
the Euclidean and the Minkowski space wave operators. They are denoted by P'^i'd) and P^(??), 
with 1? E f ]) and are defined by 



and likewise 



which can be verified by explicitly writing out the individual terms 



x^*ef{x) = ix^- i9xdt-\e^df) fix) , 

t'^0f{x) = it'+\etd,-\e^dl)fix) . 



(2.25) 




= e 




(2.26) 



P\^) = e (cos(^) (Pf)i - i sin(^) (p2)o) 



^These identities are taken in the multiplier algebra corresponding to the Schwartz space <S(R^). 
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Using (2.25) one easily checks 

P2(T?) fix) = Hii}) *2/^ fix) and P\i}) fix) = fix) *2/e ^(^) , (2-27) 

where 

Hi-d) :=E^x^+ e^'^t^) . (2.28) 

The wave operators (2.26) relate both signatures, with -i? = corresponding to Euclidean 
signature and = it ^ to hyperbolic signature. In the limit E ^ one easily verifies that this 
regularization reduces to the i e-prescription for the usual Klein-Gordon operator. Hence it can 
be regarded as a generalization of the i e-prescription to the case with an external electromagnetic 
field. To distinguish both schemes we will call this alternative prescription the -regularization. 

Using the Weyl-Wigner correspondence, the eigenvalue equations of our original operators can 
be represented on the space of Weyl symbols by 



(2.29) 



p\^) f^l-Hx) = w [h(^) f(f„-)] ix) = aL^-) f^l-Hx) , 
p'm fLl'H^) = w [f^i-) m] (^) = A(f-) fL'tHx) , 

with & = W-^i/i^"^] and the Weyl symbol 

H(^) = i (w-i[^/2i^x]2+ e2'^W-i[V2i?t]^) = ^ (p2+ e^'^'q^) . (2.30) 
The eigenvalues will turn out to depend on E and only through the combination 

Ei):=Ee'~^, (2.31) 

which explains the notation. The Hermitian symbols [v^-Ea;] = p and [-v/2i^i] = q 
obey the commutation relation of the Heisenberg algebra 

[q,p] = 2^2 w-i [t^2/EX-x^2/Et\ =4i^ , (2.32) 

where we used the fundamental property 

W[f]*2/i?W[g]=W[fg] (2.33) 

of the Weyl-Wigner correspondence. 

The operators H('!9) for -d € ( — ^,^) are known as complex harmonic oscillator Hamiltonians. 
When defined on 5(M) they have a discrete spectrum given by 



c7(H(i?)) = {A(f^) = 4^^ (m + 1) I m G Nq} . (2.34) 

The spectrum of H('!9) and its eigenoperators fmn will be investigated in §3. The simultaneous 
eigenfunctions of P'^i'd) and P^(??) are given by the Wigner transformation 

/(f„'')(a:) = W[f(i-)](^) . (2.35) 

These functions are calculated explicitly in Appendix A. They have an exponential decay for x, t — >• 
oo and are Schwartz functions. This is in contrast to the functions one obtains in the limit — t- it ^, 
which are tempered distributions and have been found in [FS09]. As shown in [Fisll, App. D], 
these functions span a dense subspace of L^(M^), thus every square- integrable function on can 
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be expanded pointwise into functions lying in the span. In addition, they fulfih the important 
projector property 



f^n'^ ix) fja'^ (x) = \ - 6^, /it ^ (x) . (2.36) 



The generahzation of the ^^-regularized wave operators to D = 2n dimensions are given by 

K^{'d) = e'^{ cos{'&) - i sin(t?) K^) and k^{'&) = e'^ { cos(i9) Rf - i sin(i?) K^) , (2.37) 

which again split up into two-dimensional wave operators defined by (2.12)-(2.14) and (2.22). In 
D = 2n dimensions the components {Pf)k and (P^)fc_i, and likewise (Pf)fc and (P^)fc_i, differ only 
by a sign for A: = 2, . . . , n up to a relabelling of the coordinates. We thus have 

n n 

K2(^) = p2(^)+ e^i'' ^(p2)^, and K2(^) = p2(^) + e^*'^ ^ (p2)^ (2.38) 

k=2 k=2 

according to (2.26). The eigenfunctions of the operators (P^)fc and (P^)fc are just Landau wave- 
functions. What remains is to find the eigenfunctions of the remaining parts of the wave operators. 
Since all of these two-dimensional differential operators commute, the total eigenfunctions will be 
a product of the individual two-dimensional eigenfunctions. 

2.3 Perturbative Quantum Field Theory 

Without loss of generality we will choose in the following always -d > Q and define = ^ — k > 
for small k > 0. Denoting 

(2.39) 

the regularized LSZ model is defined by the classical action 

s\Sz = j 'A*(^) (a(K2-^2)^ + (l-a)(K2- A) 0(^) 

-g{a j {(f)* -kQ 4> -kQ 4>* -kQ (t)){x) + P j d^x {4>* kQ (j)* -kQ (f) k^ 4>){x)^ , 
and the regularized Grosse-Wulkenhaar model by 

- g j dPx {4> kQ(/)kQ(j) kcQ (j)) (x) . 



(2.40) 



(2.41) 



The Minkowski space duality covariant noncommutative quantum field theory of the regularized 
LSZ model is defined by the partition function, which is the generating functional obtained by 
adding external sources J{x) and J*{x) to the action (2.40) with 

Z[J,J*]= lim^JV j V<j)V(t)* exp(^icSts^ + j x J* {x) (t){x) + j x (jf {x) J {x)^ (2.42) 

and analogously for the real Grosse-Wulkenhaar model, where is a normalization constant. The 
precise definition of the path integral measure is not required to determine Z[J, J*] perturbatively, 
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since only the vanishing of the integrand for |(^| — t- oo is needed to find a functional differential 
equation for the partition function via formal integration by parts in field space; this is ensured by 
the t?-regularization. The "free" partition function Zq\J^ J*] := Z[J, J*]|^^q is then the solution of 

^hm^ [a - ^?), + (1 - a) (K^^ - /x^).) ^^f,!^^^ = i J{x) Z,[J, J*] , 



lim^ {a (K^ - + (1 - a) (K^ - ^,\) ^^f/^J^ = i Z,[J, J*] 



(2.43) 



given by 



^o[J,Jl = jhn exp(i j d^x J d^y J*(^) A(-'-)(a;, y) J(y)) , (2.44) 

with A^'^''^\x,y) the regularized dressed propagator defined through the equation 

{a (K2 - + il-a) (K2 - ;,2)^) A(«-)(^, y) = 6{x - y) . (2.45) 

An explicit expression for /\^''^''\x,y) will be derived in §6.1. The full interacting quantum field 
theory is given by the partition function 

Z[J,J*] = Jnn AAexp[i5int(^,^)]exp(i j d^'x j d^y J*(a;) A(-'-)(a;, y) J(?/)) (2.46) 

leading to a perturbative expansion in Feynman diagrams corresponding to the interaction part 
5int[(/', of the action (2.5) and the dressed propagator ll^'^''^\x,y). The corresponding Green's 
functions contain products of distributions and have to be regularized; this is described in §5. For 
real scalar fields we get 

Z[J] = \im^ AAexp [i5int(^)] exp ( ^ j d^a^ j d^y J(a;) AW(aj, y) J(y)) , (2.47) 

with regularized dressed propagator ll^'^\x,y) given by (2.45) for a = ^. 

In the following we will construct dynamical matrix models representing the regularized duality 
covariant quantum field theories. In fact, our construction has killed two birds with one stone. 
First of all, we have regularized the wave operator such that no zero eigenvalues occur, so we can 
invert it to get a unique propagator. On the other hand, we have also found a discrete spectrum for 
the regulated wave operator which, together with the projector relation (2.36), is needed to define 
a proper matrix model formulation of the quantum field theory. This is in marked contrast to the 
usual i e-prescription that gives a regulated wave operator D^^^ = cr + (1 — cj) — + i e, which 
simply amounts to adding the constant i e to the continuous spectrum of the electric part of the 
wave operators, but otherwise leaves its continuous character unaltered. A perturbative quantum 
field theory amenable for the continuous basis approach with functions Xpq is analysed in this way 
in [Zahlla]. 

In the following we shall address the following questions: 



• What is the interpretation of the sign of i9? 

• Which propagator do we obtain in the limit — )• it ^? 

• Is it possible to prove duality covariance for i? 7^ it ^ at the quantum level? 

• Are the Feynman diagrams finite in the limit — )• it ^? 
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In §4 we shall argue that, like the ie prescription at ii^ = 0, the i?-regularization is related to 
causality, and flipping the sign of the regulator corresponds to interchanging the Feynman (causal) 
propagator with the Dyson (anti-causal) propagator, i.e. interchanging the particle and anti-particle 
descriptions in a background electric field. The proof of duality covariance at the classical level 
follows easily from the Euclidean and Minkowski space proofs given in [LS02] and [FS09], respec- 
tively. The spacetime metric plays no role in the proof, which relies solely on Fourier expansion 
techniques and Gaussian integrations; in §5 we will calculate the partition function of duality co- 
variant noncommutative quantum field theories, and hand in the proof of the duality invariance at 
quantum level. 

3 Dynamical Matrix Models 

In this section we will work out the matrix model representations of the perturbative quantum 
field theories defined above. In §2 we used the Weyl-Wigner transformation to map the eigenvalue 
problem for the ^^-regularized wave operators to that of the complex harmonic oscillator. Below we 
investigate its spectrum and eigenfunctions, and construct the appropriate generalizations of the 
Landau wavefunctions. Using their Fock space representation, we will finally arrive at the matrix 
model representation for the two-dimensional classical models and their corresponding quantum 
field theories. The generalization to higher dimensions is also presented. 

3.1 Complex Harmonic Oscillator Wavefunctions 

We will begin by investigating the spectrum of the complex harmonic oscillator Hamiltonian H('i?) 
defined in (2.30), with commutation relation (2.32) and positive real frequency E £ ]R_|_, which 
turns out to have a discrete spectrum (2.34) resembling the usual harmonic oscillator spectrum 
rotated into the complex plane by a phase factor e . Since q = [^/2Et], it is natural to 
work in the representation defined by the eigenbasis of q such that 

{q'\q\q) =V2Eq{q'\q) and {q \ P \q) = - i Vs dg {q' \q) (3.1) 

and thus 

{q'\Hi^)\q) = {-4d'g+Ejq'){q'\q) (3.2) 
with the condensed notation (2.31). Firstly, note that the eigenvalue differential equation 

(-40,2 + Ej q') f^-\q) = AE^ (m + 1) /(f-)(g) (3.3) 

is fulfilled for complex values if fm ^ (q) represent the usual Hermite oscillator wavefunctions 
fm (q) with the complex frequency substituted for E^ 

fL'^H^) = (:^;^r^)''' e-^-V4^^(y^,) , (3.4) 
V 2™- m! \'2-K ' 

2 2 

where Hm{z) = (—1)"* d^e~^ are the Hermite polynomials. These functions will be called 
complex harmonic oscillator wavefunctions, as a generalization of the harmonic oscillator wavefunc- 
tions to complex frequencies E^. They possess an exponential decay due to the Gaussian factor, 
and are thus Schwartz functions on M for I"!?! < ^. We expect that by continuity, for \-d\ small 
enough the eigenvalues of the complex harmonic oscillator Hamiltonian are given by the set (2.34); 
the values (2.34) are indeed the eigenvalues of H(i?) for < ^ [Dav99]. 
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The complex harmonic oscillator wavefunctions (3.4) are not orthogonal, and thus do not serve 
as a usual Hilbert space basis for 5(M). But together with their complex conjugated functions 
and for lHe(£^^) > 0, they constitute a bi-orthogonal system with respect to the L^-inner product 
(— |— ). This means that the two sets of functions {fm^^)^^^^ and {fm'"^^)^^^^ with nonzero 
and D\e{E^) > fulfin 

</f -"Vi'''^> = 7^9 /f ^^(g) f^'Hq) = Snm , (3.5) 

which follows immediately from the orthogonality of the Hermite functions on M by a deformation 
of the integration contour to a straight line from — oo e to +oo e This rotation is possible due 
to the Gaussian factor in the integrand, ensuring an exponential decay for 9^e(£'^) > 0. 

In addition their linear span is dense in L^(M), which means that every square- integrable function 
on M can be approximated pointwise by a linear combination of these functions; the proof can be 
found in [Fisll, App. D]. But the series which occur are not convergent in the L^-norm and thus 
do not build a Riesz basis [Dav99, DK04]. 

To ensure the applicability of this basis to arbitrary quantum field theories, however, one has 
also to be able to deal with scalar products and (tempered) distributions. The problem of uniform 
convergence for < ^ might be circumvented by considering a smaller space than the Schwartz 
space, like the space of smooth functions with compact support, or by considering the Sturm- 
Liouville problem for the complex harmonic oscillator Hamiltonian on a finite interval [—L, L] in M; 
the expansion on 5(M) might then be defined in some limiting procedure. The applicability of the 
Gel'fand-Shilov space of type 5°(M) C 5(]R) with a = ^ as an appropriate dense subspace of fields 
on M is discussed in [Fisll, App. C], see also [Sol07b, Sol07a]; but the question of which precise 
spaces of functions this complex oscillator basis is applicable in a distributional sense is still an 
open problem. 

Although it would be desirable to have a general rule which tells us for which functions the 
matrix basis is applicable, for a given field theory it suffices to derive the asymptotics of the matrix 
space representation of the corresponding propagator in order to ensure the convergence of the sums 
in Feynman diagrams; this is investigated in §6, but at present it is an open issue. In the following 
we will use the matrix basis to derive the propagators of the various field theories and find that they 
coincide with the position space propagators in all cases for which results are already known in the 
literature. In Appendix B the one-loop effective action of the Klein-Gordon theory in a constant 
background electric field is calculated with the help of the matrix basis, and shown to also coincide 
with the known results. By picking up the regularization scheme imposed on the position space 
propagator in the Euclidean case, which effectively cuts off the matrix element summations at some 
finite rank N, the occuring Feynman diagrams of the ??-regularized field theories are well-defined 
and duality covariant. Whether or not new divergences arise in the limit N ^ oo remains to be 
investigated. 



3.2 Complex Landau Wavefunctions 

We will now construct the complex Landau wavefunctions fmn \ defined by (2.35), through Wigner 
distribution of the tensor product of two complex oscillator wavefunctions fm^^ ■ We will also derive 
a "ladder operator" type construction, which allows us to obtain the matrix model representation 
of the duality covariant field theories. For the moment we set 9 = 2/E. 

We will first relate the ordinary and complex harmonic oscillator wavefunctions using complex 
scaling methods. Introducing the Hermitian scaling operator 

V(^)=exp(- A(pq + qp)) (3.6) 
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we see that 

\/{^)q\/{^)-^ = e'^/^q and V(^?) p = e " ^''/^ ^ . (3.7) 

The complex harmonic oscillator Hamiltonian (2.30) is thus related to the ordinary oscillator Hamil- 
tonian by 

H(^?)= e'''\/(^)l(p2 + q2)V(^)-i , (3.8) 

while the complex eigenfunctions can now easily be obtained from the orthonormal oscillator num- 
ber basis states \m), m G No, where (m|n) = 5mn and {q\m) = fm (<?) are the ordinary harmonic 
oscillator wavefunctions 

fL'H^) = e-^"'/^ h^{^,) . (3.9) 



2™ m\ V2tt 
By noting that 

H(i9) \/{'d)\m) = e '^\/{^) H(0) |m) = e 4^ (m + i) V(??) |m) , (3.10) 
the corresponding eigenfunctions are related to the oscillator wavefunctions by 

fi^'HQ) = {Q\fi^'^) := {q\\fm\m) = e'''" f^r^^e'^^'q) . (3.11) 

Clearly the left/right eigenoperators of H(i?) are tensor products of the form 

= \l^\im^){n\\l{-^) , (3.12) 
and the complex Landau wavefunctions are thus given by 

f^n'\^) = [\/(^)|m)(n|V(-^)] (x) . (3.13) 

The normalization has been chosen such that 

J d^x f'i^\x) = J dt J dx J dk e'^^'"'^ {t + k/2\\l{'d)\m){n\\l{-^)\t-k/2) 



where we have used the explicit representation for the Wigner transformation (2.17). From (2.17) 
we also see that complex conjugation yields 

f'^n^\xr = ^ I dk e-''"'-/^{t + k/2\\/{-^)\n){m\\/{^)\t-k/2) = fii-'>\x) (3.15) 
and the projector property takes the form 



'r^*2/Efif'^){^) = ^\N\\/mm) {n\k) (/|V(-^)1 {X) = J ^ . (3.16) 



mn '^2/EJkl 



E 




Together with the normalization condition this implies the bi-orthogonality of the system of complex 
Landau wavefunctions with respect to the L^-inner product 

ifLlMf''^) = I d^x fii-^\x) fif-^\x) 

= I d'x {Al-'^*2/Efif~'^){^) 

= j d^X 5mk fl!l-'\x) = 5mk5nl. (3.17) 

The explicit expressions for the matrix basis functions are given by 
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Proposition 3.18. The complex Landau wavef unctions fmn\x) for m, n G No are given by 

p|m— n|/2 



f{Ei})(f , 



^ j^^min(m,n) 



E / min(m!, n\) 



vr y max(m!, n!) ^ 



,-i^.x!,-'xL^V2(^W y--n| l^nl (^^J^)^W) (3.19) 

V — sgn(m— n)/ mm(m,n) V + — y ' v / 



(m,n) ' 



where 



t± ie-'^^x 



(3.20) 



are complex light-cone coordinates and L!^{z) are the associated Laguerre polynomials. 



The proof of Proposition 3.18 is found in Appendix A. Setting = this result coincides with 
the weh-known expression for the Landau wavefunctions in the Euchdean case. The coordinates 
(3.20) continuously interpolate between the complex coordinates tit ix in the Euclidean case ■!? = 
and the light-cone coordinates tzizx'm the hyperbolic case ?? = zt ^ . Since 



E,') x \ x _ 



^)=E[ cos(??) (t^ +x^)+i sin(i?) {t^ - x^)) 



(3.21) 



we see that similarly to fm^'^ these functions are Schwartz functions only for \'d\ < ^; in particular 
they belong to the Gel'fand-Shilov spaces 5"(M^) for all a > ^. At = it ^ they have a polynomial 
increase and are thus tempered distributions. 

The Fock space representation of the harmonic oscillator wavefunctions has a counterpart in the 
complex scaled version, which will prove very useful in the explicit determination of the dynamical 
matrix models. For this, we note that 



V(??) 

1 



(at) 



|0)(0| 



V(i?)-^ 



(V(^)atV(^)-i)'"|/f-^)(/o^^-^| (V(^)aV(^)-i) 
\/m\ n! 



(3.22) 



with a 



7= (q + ip). We can use the relations (3.7) to get 



V(??)atV(??)-i 
\/(i9)aV(??)-^ 



1 



/8E 
1 



ie-^^/2^) 



ei''/2q+ie-i^/2p) . 



Since W^^ [V2Et] = q and W^^ [^2 E x] = p we find 



W 



V(??)atV(t?)-i 




and 



W 




4 ^+ ' 



(3.23) 



(3.24) 



where we used the complex light-cone coordinates (3.20). The corresponding derivatives are given 
by 



5iT ie'''a. 



(3.25) 



with x^^ = 2 and d^^x^^ = 0. The matrix basis functions on M? can now be obtained via the 
Weyl-Wigner correspondence 



J mn \ 




An 
1 



W 



(3.26) 



V ml nl 



4 



''2/E 1 




w 



4 
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We define "ladder operators" through the equations'^ 

^ xL^^ ) -k2iE g{x) = a^E^^g{x) and ) 9{x) = a~j^^^g{x) , 



^(a;) *2/£; ( Y ^ ) = 6+j^)C/(a;) and g{x)-k2iE (^^^ x''^^^ =b^^^^g{x) 



(3.27) 



The differential operators on the right-hand sides of these equations can most easily be obtained by 
expressing the star-product in terms of the complex hght-cone coordinates. Inverting the relations 
(3.20) and (3.25), after a bit of algebra we find that the ladder operators are then given by 

f.„4(v/^4"'.\/i4"') and ,J=,^,4(^.f,y^af 



"rE,) = t(v/^*4"'=F\/i^"-' '^fE., = -A■JE,:^T^^IiraT]■ (3.28) 

and that they fulfill the commutation relations 

['\e,) ' = 1 and [h-^^^ , 6+ = 1 , (3.29) 

with all other commutators equal to zero. As expected, we arrive at the usual Euclidean case from 
[LSZ04] when substituting by E. 

The ground state wavefunction is determined by the differential equations 



^-iE,)fm'\^) = blE,)f!>o'\^) = (3.30) 



plus the normalization condition (3.14) with m = n = 0, which has the solution 



( F ) 

The wavefunctions fmn have the ladder operator representation 



(a+ (b+ V 

/(f„-)(^) = f^^^'\x) . (3.32) 



It immediately follows that 

^iE,)f^nH^) = V^f!n-U^) and aJ,^)/(f/)(^) = V^^fi,f:Ux 

b^Ej^n^K^) = V^fSnU^) and b+^^^f(^^^\x) = V^fS:U^) 



(3.33) 



We will use these relations to obtain the desired matrix model representations. Note that, by [FS09, 
Lem. 5], the problem of the right test function space is the same as in the complex oscillator case 
of §3.1; we can relate the subspaces of Gel'fand-Shilov spaces 5^(M) to subspaces of 5^(M^) via 
Wigner transformation. 

3.3 Matrix Space Representation 

Using the Fock space representation of §3.2, we will now derive the matrix model representation of 
the classical regularized actions of the LSZ model. In the following we denote 

/™n := f!ni'-'^ (3.34) 



^Since (a^^^)^ / "^(Bs) ^^'^ ^''CEii)^^ ^ ^(Eh) ^^^y ^''^ strictly speaking not ladder operators in the conventional 
sense, but we will nevertheless refer to them as such. 
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with "!? = ^ — K and va. general; in this case the complex Landau wavefunctions diagonalize 

the interaction part of the action, but not necessarily the free part of the action. 

We expand the scalar fields in terms of the complex Landau basis^ 

oo oo 

0(a^)= /mn(a;)Cn and (t^ixf = ^ /^n(a^)Cn, (3.35) 

m,n=0 m,n=0 

where the complex expansion coefficients are given by 



Cn = (/mnl0) = j ^ ^ f^^.^{x) ^{x) and Cn = (/mnl</'* ) = J d'x (P{x)* 

(3.36) 

with (pmn — i'^nm)* ■ The free parts of the actions can be deduced from 

Lemma 3.37. The i!) -regularized wave operator of the 1+1 -dimensional LSZ model in matrix space 
is given by 

dI::1 = (-e->2^2iii±^(m + n + l) + li^(n-m))5™,5„, 

■)2 



+ 2i — - — \^^/nm6m,l+lSn,k+l + + 1) (m + 1) (5^,/-! 5n,fe-i j (3.38) 
with frequencies = E 9/2 and = {2a — 1) 0. 
Proof. The wave operator is defined in the matrix basis by 

Dl:nt = I d'^/mn(^) (ae-p2(f-K) + (l-a)e-p2(f-K)-e--/i2) /^,(^) . (3.39) 



One has 

29 

+ {9^ - 4) (a+/,^) 6+/,^) + a-/,^) ft"/,^))) , (3.40) 

together with a similar expression for P^(t?) with ci^2/e^) ^(z/e^) interchanged. These formulas 
are verified directly by inserting (3.28). The matrix space representation of the partial differential 
operators P^('!?) and P'^{i!}) away from the self-dual point can be obtained from (3.40) with the help 
of (3.33) to get 

Pln;M = ^{{2 + E9f{m+'^)6„uSnk + {2-E9f{n + ^)6rni6nk (3.41) 



+ [9"^ E"^ - 4) {y/nm5m,i+i 5n,fc+i + \/(n + l)(m + l) 5m,i^i 5„,,fc_i) 



and 



PL;hW = -^[{2 + E9f{n + \)5mi5nk + {2-E9f{m + \)5mi5nk (3.42) 



+ (02 ^2 _ 4^ (^ATm (5^ 5n,k+i + V(n + l)(m + l) 
which can be combined to give (3.38). □ 



This expansion is defined for L -functions in a limiting procedure which can be found in [Fisll, App. D] 
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For the LSZ interaction terms, we use the projector property (3.16) to get 

/mini '^0 fm2n2 ^9 fm^n^ '^0 fm^rn (27r 0)'^/^ ^n2mz ^risnn f^,n, ■ (3.43) 

The regularized LSZ model in two-dimensional Minkowski space can then be represented in the 
matrix basis as 

S\Sz= E Cn^LXi'^rfc-^ E («CnCfc^feZ'/'L+/?Cn^nfc'A^/'^L) • (3-44) 
m,n,k,l m,n,k,l 

As Si oiie-nicitrix model with infinite complex matrices (pi^ — 

{4'mn) m nsNo representation reads 
Slsz = ^ '^^((2 - ^ + 8aeE) f 0, + ((2 + 9E) -8a9E)^^£ 4>l, 

+ i {9^ - 4) {<pl, rt 0, r + rt T)-2ee~'^^? </.^ (3.45) 
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27r6' 

with the diagonal matrix 

£mn = 4i (m + i) 5ran (3-46) 

and the infinite shift matrix 

Tmn = \/n-l Sm,n-1 ■ (3.47) 

Using the perturbative solution (2.46), the duality covariant field theory can be defined pertur- 
batively in the matrix basis by the partition function 

Z[J,J] = hm AAexp(-^ J] f ) (3.48) 

with , JJ^„ external sources in the matrix basis and the propagator A^^°^].; defined as the inverse 

of n^''''^) 

E ^rnn-]cl ^Zfc;'sr ~ E ^nm-Jk ^kl;rl ~ ^ns ■ (3.49) 

An explicit expression for the propagator ^^^^-ki ^'^^ derived in §6.2. The modified Feynman 
rules are presented in the double line formalism and are exactly as in the Euclidean case [Sza03, 
GW05b]. The generically non-local propagators are represented by double lines with orientation 
pointing from 0* to as 

n ^ k _ aCk.'^) 

> nmdk ' 

m I 

The two interaction terms 0* 0*5i 0* -kg and 0* *g 0* ★g 0*e are represented by different diagonal 
vertices given respectively by 
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TT 

2w 6 ^"rnp ^nq ^kr ^Is and — X---^ r--^ — — 2^ 9 ^nq ^kr ^Is 



Restricting to one of these interactions reduces the number of possible diagrams for the complex 
matrix model. 

For real fields, one can apply Lemma 3.37 by setting cj = ^ to immediately get 

Lemma 3.50. The •& -regularized Grosse-Wulkenhaar wave operator in 1 + 1 dimensions has the 
matrix space representation given by 



1^2-1 



+ 2i — - — yVrim6m,i+iSn,k+i + \/(?^ + 1) (m + 1) Sm,i~i 6n,k~ij (3.51) 
with frequency J7 = E9/2. 

The Minkowski space Grosse-Wulkenhaar action in the matrix basis then reads 

4w= E (^Cn^2;H'A^/- ACnCfc'/'fci^L) (3.52) 



m,n,k,l 

with (f)'^^ = 4>nm'^ it thus corresponds to a Hermitian one- matrix model. From the perturbative 
solution (2.47) the partition function in matrix space is given by 



m,n,k,l 



m,n,k,l m,n,k,l 



where the propagator A^^.^^ is the inverse of Dl^n ki ^^"^ represented by the unoriented double 
line 

» nm\lk 

m I 

The vertex of the (j)*^ interaction is given by the graph 



2-k9 ^"m-P ^nq ^kr ^Is 



TA 



Since the vertex is unoriented there are as many diagrams as in the LSZ model with both parameters 
a and /3 turned on. 
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3.4 Generalization to Higher Dimensions 

The spectra of both partial differential operators in (2.38) in generic D = 2n dimensions are given 
by the set 

4£;e^''(/o + i) + E 45fe e2^''(Zfc + i) ^o, ^i, • • • , ^n-i G No L (3.54) 
k=i ^ 

where the eigenfunctions are products of the complex Landau wavefunctions fmono from §3.2 and 
the ordinary Landau wavefunctions fm^n^ , so that 

n-l 
k=l 

with Xk = {x'^'' , x'^'^'^^) G M? for k = 0,1,..., n-l, x = {x'^) = (x^, a;^ . . . , x'^) G M^, m = 
{mk),n = {uk) G N[}, and = (Ei), Bi, . . . , Bn-i) G C+xM^j:"^ where C+ := {z G C | me{z) > 0}. 
The star-product of two multi-dimensional complex Landau wavefunctions (3.55) with respect to 
the deformation matrix (2.9) decomposes into star-products of Landau wavefunctions depending 
on Xk for A; = 0, 1, . . . , n — 1. If in addition E = 2/6q and = 2/9^ for A; = 1, . . . , n — 1, then 

{f^ri^ *e /if (^) = (3.56) 

with 5mn = Yik ^ra^nk ■ 

The generalization of the matrix model representation to higher spacetime dimensions is now 
straightforward. To conform with our previous conventions we again set = ^ — k > and use 
the notation 

n-l 

fmni^) = /mono (^^o) JJ fminl\^k) ■ (3.57) 

k=l 

The functions f!^^ are arranged so as to diagonalize the interaction part but not necessarily the 
free part of the action. The scalar fields on M.^ are expanded in the complex Landau basis 

H^)= E /mr^(^)'/'^n and <A(^)* = J] f^J^)^^mr. , (3-58) 

where the complex expansion coefficients are given by 

if^nU^) = J d^'x f^^{x) <l^{x) and -^-^^ = {f^^JcP* )= I d^x f^^{x) c/^ix)* . 

(3.59) 

The matrix space representation of the LSZ model in D = 2n dimensions away from the self-dual 
point can be obtained by comparing the operators (2.38) with their two-dimensional constituents, 
and their matrix representations given by (3.41)-(3.42) together with their Euclidean counterparts 
for = [LSZ04]. The matrix space LSZ wave operator is thus the sum of the two-dimensional 
Minkowski space wave operator given by (3.38) plus n — 1 copies of the massless Euclidean wave 
operator for -d = [LSZ04] times the factor — e~'^. Noting that the massless LSZ wave operators 
in Euclidean and Minkowski space differ only by a factor of the imaginary unit i , we can write 

n-l 

= l^mono;fco«o ~ 1^ ^m^n^-^k " ^ ^^ ^rr^l Sr^k , (3.60) 

i=l 



■^mn \Jmnl 
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where m = {mo, mi, . . .,mn),n = (no,ni, . . . ,nn),k = {ko,ki, . . . ,kn),l = {lo,h, ■ ■ ■ ,ln) G Ng and 
'^mrvki two-dimensional massless Euclidean LSZ matrix space wave operators 

'^Jliki = (2^^(m + n + l) + ^(n-m))5™z5„fc 



+ 2 — [Vrim5m,i+i Sn,k+i + \/(n + 1) (m + 1) Sm,i~i 5n,fc-i I , (3.61) 



with frequencies = E9q/2 = Bi9i/2 and = {2a — 1)0.. The 2?i-dimensional regularized LSZ 
action is then given in the usual matrix space form 



= E <t^^m..DtSuA (3.62) 

m,n,k,l 

Vdet(27re) ^^^^ 

Every other result of this section (and ensuing ones) can now formally be generalized to higher 
dimensions by substituting multi-indices m, n, . . . G Ng for the usual matrix indices m,n, . . . G Nq. 



4 Causality 

In this section we will treat the problem of determining the causal propagator of the duality 
covariant field theories in Minkowski space. Problematic for this issue is the lack of time translation 
invariance, which allows for transitions that violate energy conservation; this manifests itself in an 
instability of the vacuum with respect to production of particle-antiparticle pairs. We review how 
the standard techniques must be altered to take care of these features. We will then examine the 
corresponding propagators which one obtains by removing the iJ-regularization. 



4.1 Causal Propagators 

The way we chose the propagator of the Minkowski space theory was to find the analytically 
continued propagator of the Euclidean case. This also brought about the possibility of finding a 
matrix space representation. In the following we will show that the propagator which is prescribed 
by the iJ-regularization is the causal propagator of the duality covariant quantum field theory. For 
this, we will first review the connection between regularization and propagators by describing the 
eigenvalue representation, and the related operator extension method. 

The free partition function Zq[J, J*] of a complex scalar field theory is defined as the vacuum- 
to-vacuum amplitude 

Zo[J,J*] = (O,out|0,in)'^"^* , (4.1) 

where |il,in) and (0,out| are the vacuum states at time instances tin and tout of the quantum 
field theory defined by the free action So[(j),(j)*] in the presence of the sources J and J*. Using 
Schwinger's action principle, one can show that causality implies the relation 



6^ logZo[J,J*] 
6J*{x)5J{y) 



J=J*=0 



(O,out|r(,^(a;),^(y)t)|0,in) 

(O,out|0,in) ' ^ ' ' 



where 0(a;) is the second quantized field operator, T denotes time-ordering with respect to the time 
variables and ip , and |0, in) and (0, out| are the in- and out- vacuum states for J = J* = which 
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in the presence of further interactions are taken in the interaction picture where the field operators 
satisfy the equations of motion obtained from varying 5o[</>, <^*]. For field theories which allow for 
spontaneous particle-antiparticle pair production, like the covariant models we are considering, the 
in- and out- vacuum vectors are in general not dual to each other. Thus a non-trivial vacuum- 
to-vacuum probability |(0, out|0, in)p < 1 may occur, since |(0, out|0, in)p measures the vacuum 
persistence which is equal to 1 only if no spontaneous pair creation occurs. 

The tempered distribution defined by the right-hand side of (4.2) is known as the causal propa- 
gator and will be denoted as i Ac{x, y), where the imaginary unit has been factored out to conform 
with our previous conventions. Quite generally, for a Klein-Gordon field which may be free or 
moving in an external background which preserves vacuum stability, the expression (4.2) may be 
evaluated through the expansion 

iA,(a;,y) = T(xO-yO) [ dm{u) 4>i+\x) 4>M (^y)* + T{y^ - x'^) [ dm(z.) 4-)(^) ,^i-)(y)* (4.3) 

Jc Jc 

with r the Heaviside distribution function, {(pl^^) a complete set of distributional solutions of the 
classical field equation with positive and negative frequency, respectively, and dm(z^) a suitable 
measure on the set C of all quantum numbers u parametrizing the space of solutions. One can 
check that the distribution (4.2)-(4.3) propagates particles (positive frequency solutions) forward 
in time and anti-particles (negative frequency solutions) backward in time. This is the imprint of 
causality and lends the causal propagator its name. 

The situation is more complicated if the background field spoils vacuum persistence. A typical 
example is ordinary quantum electrodynamics in an external field which allows for pair creation. 
Crucial for the canonical quantization scheme and for the expression (4.3) is the existence of a 
complete set of classical solutions which have definite positive or negative frequency for all times. 
However, such a set of solutions only exists if we are working on a stationary spacetime, i.e. 
a spacetime which admits a global timelike Killing vector field [DeW75]. In our case, there is 
no such vector field due to the absence of time translation symmetry; production of particle- 
antiparticle pairs manifests itself in an inevitable mixing of positive and negative frequencies at the 
level of solutions to the field equations. The requisite methods in this case have been developed 
in [Git 77, FG81]. 

Since the asymptotic Hilbert spaces in the remote past and future (if they exist) are differ- 
ent, there are two complete sets of solutions denoted (<^iy^^^) and ((/)j^(_|_)), having definite posi- 
tive/negative frequency at times tin and tout, respectively, which are the equivalent of the posi- 
tive/negative frequency solutions above in the infinite future and past, i.e. in the limits tin — ^ —oo 
and tout — ^ +00. The generalization of the expansion into classical solutions (4.3) then reads 

iA,{x,y) = r(xO-y°) / dm(^) / d?n(z/) 0^(+)(^)a;+(/i|z.) ,/.,(+)(y)* 

Jc Jc 

+ r(2/0-x°) / dmifx) [ dm{u) <P,(^^){x) uj- ^^^-^yT . (4.4) 
Jc Jc 

Here uj^{ij,\u) is the relative probability for a particle/anti-particle to be scattered by the vacuum 
in the external field, given by a generalized Wick contraction of creation-annihilation operators on 
Fock space which appear in the mode expansions of the in- and out- field operators, and which create 
the in- and out- particle/anti-particle states. For a field theory with a stable vacuum state one 
has io^d^lv) = 5{ii,v) and = 4>u[±): where dm(/i) 5{fi,ij) f{^i) = f{y). This construction 
determines the propagator uniquely and is equivalent to the definition (4.2), but can be quite 
technically cumbersome to carry out explicitly; hence it is desirable to have another method at 
hand. 
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Such an equivalent method, which will prove profitable for us, is the eigenvalue representation. 
Let (p\{x) be a complete orthonormal set of eigenfunctions of the wave operator Da, of the field 
theory with eigenvalues A E a{Dx), i.e. 

Dx(px{x) = Xipx{x) (4.5) 

with 

/ d^(A) ipxix)* My) = 6{x - y) and / d^cc ^x)* ^x'{x) = 5{\, \' ) , (4.6) 

Jcr(D^) J 

where the measure d^(A) is discrete measure on the point spectrum and Lebesgue measure on the 
absolutely continuous spectrum in a{Dx) C C with J^^p ^ d^(A) 5(A,A')/(A) = /(A'). Contrary 

to the functions (p^^ above, these eigenfunctions need not solve the field equations. Decomposing 
the propagator into these eigenfunctions gives the formal expansion 

A{x,y)= [ d£{X) ifxi^r My) ■ (4.7) 



However, the potential poles at A = make this definition problematic, which reflects the 
existence of more than one propagator for a given field theory. Usually one modifies the denominator 
by adding an adiabatic cutoff A — ?■ A + ie-F(A) with small e > and a function F : cr(Da;) — ?■ M on 
the spectrum of Da;, so that 

A+ieF(A)/0 (4.8) 

for all A € (7( Da;). A Green's function for the partial differential operator Da; is finally obtained by 
taking the adiabatic limit e — )■ 0^. Equivalently, one can regularize the operator Da; — )• D^*"^ with 
lim^_i.o+ D^*'^ = Da; and solve the equation D^*'^A('^^(a;, y) = S{x — y), where limg_^Q+ A^'^\x, y) is a 
Green's function of the original wave operator Da;. 

Hence any well-defined operator which is continuously connected to the original wave operator 
and has no zero eigenvalues gives rise to a propagator for Da;. However, apart from the requirement 
of absence of zero eigenvalues of D^*"^ (or equivalently the condition (4.8)), the regularization is 
arbitrary and different regularizations may lead to different propagators. For example, in the free 
Klein-Gordon theory the choice of F{k) as a positive constant leads to the Feynman propagator, 
while F{k) = 2k^ yields the retarded propagator. In general, one cannot be sure whether one 
obtains the causal propagator unless one compares it by hand to the result obtained from (4.2). 
This is the obvious drawback of the eigenvalue method, and while the problem is easily solved in 
free field theory, it is still unsolved for the general case of an arbitrary propagator and arbitrary 
external field; in particular, the equivalence of the propagators in the different representations for 
generic electromagnetic backgrounds is still an open question. 

For the LSZ model we already encountered the two regularized wave operators D^^ given by the 
?9-regularization and the i e-prescription. The question of which propagator they lead to in the limit 
e — 7- 0"*" has been answered for the i e-prescription for several related models. For the Klein-Gordon 
field moving in crossed or parallel uniform electric and magnetic fields, or in an electric field with an 
additional plane wave background, this method gives the causal propagator [Rit70, Rit78, BFS85]. 
Since an additional uniform magnetic background should not change the pole structure of the 
propagator, the i e-prescription should also give the causal propagator in the background of a pure 
electric field. In §4.2 we will confirm that the "iJ-regularization also gives the causal propagator in 
the background of a uniform electric field along one direction. 
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4.2 Causal Propagator in Matrix Space 

Using the "sum over solutions method" (4.4), the causal propagator for a massive complex scalar 
field of charge e in four dimensions in the background of a constant electric field along one space 
direction has been calculated in [FGS91, eq. (6.2.40)] with the result 

a E _ : „ ™,. . IT,,,, /9 /"°° ds 1 



s sinhfs eE) , , ^, 

„ 1,2 ("^-9) 

/ 2 1 n II ii2 ^ / \\x± — y±\\p 
X exp — IS fi — 7 e £/ kcii — yii Ui cothfs eE) + i 

Here we defined x = (a;|| , x±) G with x± G denoting the two space components perpendicular 
to the electric field and 

x^\-Ey\\:=Ex^^e,,,yl^ , (4.10) 

where e^u is the rank two antisymmetric tensor with eoi = 1 and E > the electric field strength. 
Below we will start with this four-dimensional wave operator, with the electric part regularized as 
in (2.45), and calculate its (unique) propagator. For k — )• 0"*" we find exact agreement with (4.9) 
confirming that this is the causal propagator. The calculations performed here using the matrix 
basis are comparably simple, so that the matrix basis can be alternatively regarded as a powerful 
computational tool in ordinary field theory. 



We begin with some notation and a preliminary result. We define the symmetric bilinear form 

(-, -)^ : ]R2 ^ M2 ^ C for 1? G [-f , f ] by 

{x,y)^ = cos(t9) {x, y)E + i sin(^) {x, y)M = i e (xf y^^'' + xL^^ ) , (4.11) 

where (— ,— )e is the two-dimensional Euclidean scalar product and (— ,— )m the two-dimensional 
hyperbolic inner product. We also define the map || — ||^ : — t- C by 

||£c||| = {x,x)^ = cos(i9) i sin(^9) = e'^x'-^^x^^^ (4.12) 

with II — ||e the two-dimensional Euclidean norm (2.2) and || — ||m the two-dimensional hyper- 
bolic norm (2.1). For arbitrary two-dimensional vectors x,y £ we denote as above x ■ Ey = 
E x^^ e^u . In Appendix C we prove 

Lemma 4.13. For any a;, y G and a G C*, one has 

Y.f^r^n\^)fk^\y)oJ'' = -^e^v{-i^\\x-y\\l + {a-l)E{x,y)^-aix-Ey) (4.14) 

n=0 

X L^^[E \\x - y\\l - a(l - a~^fE{x,y)i) + {a - a^^)\x ■ Ey) . 



Now we determine the propagator of the Klein-Gordon field in four dimensions coupled to 
a constant electric field, where the wave operator parallel to the electric field is given by the 
two-dimensional ^^-regularized operator (P^ — /i'^)^. The coordinate vector is again written as 
x = (a;||,a;_L) G M^, with x± G ker[E) the components perpendicular to the electric field, and 
analogously for the momenta p = {p\\,p±) G (M^)* and derivatives = (i9||,i9_l). 

Proposition 4.15. The propagator of the ^-regularized wave operator (/^ — ^^)^ = {P^ — /^^)^ + 
{id±)'^ coincides in the limit k — )• 0"*" with the causal propagator (4.9). 
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Proof. The inverse of (K^ — ^^)^ is given by 

A(«'^)(a;,2/) = (^1— , (4.16) 

where (K^ — /x^)^ = e P^(f — i^) — e'^'/i^ + (i9±)^ with «; > has the eigenvalue equation 
(K^' - /^') e-'^-^-) = (4ii? (m + i) + llp^lll - e-' /i^^^^ 

(4.17) 

with ■(? = ^ — K. We write fif. := e~^'^ fu? for brevity, remembering that it has a small negative 
imaginary part. Using the identity 



1 

a 



oo 

IS a 



i / ds e'"" for Om(a) > , (4.18) 



we obtain 



POO p j2 

Jo J i^^y 

X g-i^/^K g-4s-E{m+i) g is||px|||-i(:cx-y±)-P± 

The sum over n is given by Lemma 4.13 with a = 1, and the resulting sum over m follows from the 
identity [Han75, eq. (48.4.1)] 

°° 1/7/ ^1/2 I +-1/2 s 

m=0 

which yields 

^ . /.oo exp ( - is/x2 - - yiilll coth(2s£;) 



27r 7o sinh(2s £;) 

X y ^^-^ g is||P±ll|-i(a'x-y±)-Px _ (4-21) 

The integration over the perpendicular momenta can now be performed by using 

^ isp^-i{x-y)p _ ^ I i^ p i(x-s/)V4s 

to get 



dp e = \l— e ^^^-w , (4.22) 



A("'i)(a;,y) = e-^-ii-'^^'ii / — ^ r (4-23) 

^ Svr^ Jo s sinh(2sS) ^ ^ 

X exp y — is i-Q — -E \\x\\ - y\\ \\^ coth(2s E) + \ - 



4s 

Taking the limit k — )• O"'', thus — ?■ ^, and substituting E ^ eE/2 to conform to the conventions 
of [FGS91], this result is identical to (4.9). □ 
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The eigenfunctions for the full regularized wave operator (K^ — /^^)^ factorize into components 
perpendicular to the electric field times the eigenfunctions of (P^ — Since the eigenvalues of 

the perpendicular momentum operators do not produce new pole singularities, we can neglect them 
in this calculation and also in the calculation leading to (4.9). This result thus easily carries over 
to the two-dimensional case confirming that the "iJ-regularization imposes causality of the critical 
LSZ model. We expect that the i?-regularization also leads to the causal propagators for a ^ 1. 

The Schwinger parameter s > introduced in (4.18) only allows for the regularizations i9 > 
and jj? — ie because of the requirement Jm{a) > 0, where the latter regularization is normally 
associated to the Feynman boundary condition on the propagator. The other choices < and 
/i^ + i e can be applied by using 

- = i / ds e ^'"^ for 3m{a) < . (4.24) 



a 

The regularization /-i^ + ie is known as the Dyson boundary condition, which leads to an anti- 
causal propagator where anti-particles travel forward and particles backward in time. This strongly 
suggests that the regularization < leads to the Dyson propagator. 

The regularization of the mass /i^ = e~"^/i^ is actually irrelevant for the analysis above. Its 
only function is to provide a continuous interpolation between the hyperbolic and Euclidean wave 
operators with the help of the parameter alone, without the need to keep track of additional minus 
signs in front of the mass term. This means that the interpretation in terms of Feynman/Dyson 
propagators for the cases ■(?—)■ it ^ still holds by regularizing only the operator P^. 



5 Quantum Duality 

In this section we will treat the problem of implementing duality covariance at quantum level for 
our field theories on Minkowski space. The i?-regularization allows us to regularize the covariant 
field theories such that the duality is preserved at quantum level. This is done in the same spirit 
as in [LS02, FS09], with the i9-regularization now being the only new ingredient. In the following 
this will be demonstrated for the two-dimensional Grosse-Wulkenhaar model. The more general 
case of the LSZ model is treated in exactly the same way. 

We only need to address how the ??-regularization affects the behaviour of quantities under the 
duality transformation. The regularized propagator with = f — k > reads 



A W^T lA - /rlfi P2 + 1 p2 _ /,2^-l|„,\ _ fLn\x) firr?\y) , . 



In Appendix D we show that the Fourier transformation of the matrix basis functions is given by 

( — i \'m—n 

nf'^n'^ik) = Ali^'Hk) = ^-^^ fLl^\E^'k) (5.2) 

with E-^k = -E-^ {k^,k°).^ Since 

^[(P2(^) + p2(^?))/(f„'')] (k) = AE^{m + n + 1) ^[/(f/)] (k) , (5.3) 

^ There is a subtle difference here between the Euclidean and hyperbolic cases. Contrary to the ordinary Landau 
wavefunctions in Euclidean space, the (unsealed) Fourier transforms of the complex Landau wavefunctions have 
interchanged indices m •<-> n and a reflected regularization parameter i9 —iS. The interchange is equivalent to time 
reversal (or parity), see Appendix D. The reflection corresponds to charge conjugation, i.e. exchange of particles with 
anti-particles; this follows from the results of §4.2, where the regularization i9 > is identified with the Feynman 
boundary condition and i9 < with the Dyson boundary condition. The specific rescalings of momenta from (R'^)* 
to in both cases, which are formally identical but differ by the signature of the metric applied, compensates for 
this difference. 
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we find that Fourier transformation relates the propagator in position space to the momentum 
space propagator even in the regularized case as 

{T T) [A^'^)] (fc, P) = ^ {E-^k , E-^p) . (5.4) 

This relation just reflects the classical duality covariance for g = 0. 

Analogously to the Euclidean case [LS02], the UV/IR-symmetric regularization now amounts 
to cutting off the matrix element sums at some finite rank by modifying the regularized position 
space propagator to 

A^^\x, y) = {x\ (i Pi + i P2 - L( A"2 I p2(^) + p2(^) 1)1^), (5.5) 

where A G 1R+ is a cutoff parameter, and L : — t- [0, 1] a smooth cutoff function which is 
monotonically decreasing with L{z) = 1 for z < 1 and L[z) = for z > 2. We adjust the matrix 
basis functions so as to diagonalize the regulated Grosse-Wulkenhaar propagator 

—^L{Ak-^E{m + n + l)). (5.6) 



2iS(m + n + l) - e-i'^/i^ 



The interaction vertices in the matrix space representation are now quite complicated; they are 
proportional to 

v''{mi,ni- ...■,mi,n^)= j d^x (/;^^„^ -kg -kg f^^^^^^ -kg /;;^„J (x) (5.7) 

with 9 7^ 2/E in general. Since for k > the complex Landau wavefunctions are elements of 
the Gel'fand-Shilov spaces 5°(M^) with a > ^, which are closed under multiplication of functions 
with the star-product, the interaction vertex (5.7) is well-defined. 

Feynman diagrams can now be obtained by taking suitable combinations of derivatives of the 
partition function (3.53) with respect to the external sources involving the regularized propagator. 
Denoting 



^AUki=--^-^''^'^^^A\^^n) , (5.8) 



they have the schematic form 



K 



ni,mi,...,nj(,m.j( k=l 

where (• • • ) denotes the contributions from products of the noncommutative interaction vertices 
(5.7) and combinatorial factors. Since the propagator A^^(m, n) is nonzero only for AE (m+n+l) < 
2A^, which at finite A is only true for a finite number of distinct values of (m, n) £ Nq, every 
Feynman amplitude is represented by a finite sum and thus constitutes well-defined duality covariant 
Green's functions in the matrix basis; this circumvents the issue of the appropriate test function 
space for the time being. By multiplying these expressions with fm^mi'^i) for i = 1, . . . , M, we get 
back the position space Green's functions with M external legs by summing over all mi,ni. They 
are also well-defined and duality covariant, since they are built from finite sums of well-defined 
covariant objects. This establishes the quantum duality in Minkowski space for the case k > 0. 
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To prove the duality covariance at k = in the same manner as above, one has to ensure that 
the interaction vertex (5.7) away from the self-dual point is well-defined. In the absence of further 
analysis, the i9-regularization should be kept while the matrix cutoff is removed, and all summations 
and integrations have been performed. Of course the limit A — )■ oo can still be ill-defined and may 
require renormalization; removing this regularization requires a good decay behaviour of the matrix 
space propagator for large values of its indices, see §6. In addition, the results from §3.1 are not 
able to exclude the possibility that even at finite k > there might be extra divergences at A — )• oo 
if we work in matrix space, stemming from the complex matrix basis itself. This, however, does not 
affect the duality covariance of the quantum field theory, which has been achieved for the Green's 
functions in position space through the regularization of the propagators in (5.5). This result is 
independent of the matrix basis. 



6 Asymptotic Analysis of Propagators 

One of the most intriguing features of Euclidean duality covariant field theories is their renor- 
malizability. We will not attempt to prove here the renormalizability of their Minkowski space 
counterparts, but start this program by deriving their propagators in position and matrix space 
representations, and studying their asymptotics. We begin by extending the formulas given in 
[GRVT06] to the hyperbolic setting, giving the propagators for the general LSZ models in generic 
D = 2n spacetime dimensions in the position and matrix bases. 



6.1 Position Space Representation 

In the notations of §3.4 and §4.2, the main result from which all causal propagators in Minkowski 
space and their Euclidean counterparts can be derived is 

Proposition 6.1. The propagator of the regularized LSZ model in D = 2n spacetime dimensions 
is given by 

^(Ka), ^ . ^ P , e-^'^' / sinh(2s^_^) . 

A^''''^^(a;,y) = - 1 e — / ds — — — r exp -j —^ixo-Eyo 

^ 2-K Jq sinh 2s^_i,) ^\ sinh(2s£^_^) ^ 



( 1 , , ,,9 ,, ,,9\ cosh(2s E'-w) ^, 

exp - - coth 2s^_^^ ico S + Wvm) + ■ x,)o f i E{xo,yoh 
\ 2 ' sinh(2si^_^j 



X 

n-l 



X 



nBk 1 / sinh(2s5fc) . 

2iT smh(2si3fc) V smh(2sisfcj 



(6.2) 



/I N /II ii2 ii2\ cos\i(2s Bk) 

X exp - - coi\i{2s Bk) Bk [\\xk\\^ + \\yk\\^) + ^:^^Yi{2sBk) ^^^^y^i^ 

with i9 = |- K>0, /i2 = e-''^^^^ E = {2a - I) E and Bk = {2a - l)Bk. 

The proof of Proposition 6.1 is found in Appendix E. We can now read off the causal propagators 
for the four-dimensional LSZ and Crosse- Wulkenhaar models. Since (— , —)-,t/2 = '^{~i —)m and thus 

II - 11^2 = ill ~ IIm' fi'^d'^ 

Corollary 6.3. The causal propagator of the LSZ model for generic a G [0, 1] in four- dimensional 
Minkowski space is given by 



A(0''^)(ic,y) = / ds 



lEB f°° g-sM -^M-^ 



E 



(27r)2 Jq sm{2s E) smh{2s B) 

/ sin(2s^) ^ sinh(2sB) „ , , 

X exp - ixo-Eyo- ixi-Byi] 6.4) 

V sm{2s E) smh(2si3j 
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with 

A E , /II ii2 II ii2 \ cos(2s£') 9 

-4m = -- cot(2s£;) (IIsoIIm + II2/o||m) + ^T^^^^j^-^(a;o,yo)M , 

= f coth(2.ij) (ll^illl + ll^illl)- 2t£f) ^^"'^'^^^^ ■ ^^-^^ 

Corollary 6.6. The causal propagator of the four- dimensional critical LSZ model in Minkowski 
space is given by 



CO 



A(0'i)(^,y) = _ AM e- -0-^3^0- i^i-Bm / ds-— ^— -— — (6.7) 



(27r)2 7^ sm{2s E) smh{2s B) 

X exp ( i E \\xo — yoWu cot(2s E) — \\xi - yiHl coth(2s B) ) . 

Corollary 6.8. The causal propagator of the four- dimensional Grosse-Wulkenhaar model in hy- 
perbolic signature is given by 

A(^\x,y) = Hds 



(27r)2 Jq sm{2s E) smh{2s B) 

/I E 
exp (--B cot(2s^) (||a;o||M + IIi/oIIm) - ^-^^^^^^ {xo,yo)i 

exp(-^i?coth(2sS) {\\xi\\l + WviWl) + w^^^^ (a;i,m)E) . (6.9) 



The Euclidean parts of the propagators here coincide with those found in [GRVT06] after 
suitable redefinitions of parameters. 



6.2 Matrix Space Representation 

Below we set 6*0 = ^1 = • • • = 0n~i =: ^ for simplicity. 

Proposition 6.10. The matrix space propagator for the 2n- dimensional regularized LSZ model in 
Minkowski space is given by 

AS:i+.,^.. = -e-A^^d..--(--+i)-^J:r.-.^(--+iM+^ (6.11) 

n-1 



X a'-"-* ('i) TT ('i) 

^mo,mo+"o;'o+oo,'o^ ' ii mi,mi+ai;li+ai,li\'^J 



i=l 



with hyperbolic part 



mm{m,l) ^_ie'''u/'i ^_ ; g i « ^^ m+i— 2n 



'^m,m+a;l+a,li-^) = 2 .^xa+m+/+l (6.12) 

n=max(0,-Q) U " \l+n)2 ° ) 
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and Euclidean part 

min{m,l) 



^m,m+a;l+a,li^) = X] 7 \ a+m+l+1 (6.13) 



«=max(0,-a) ( 1 — (-^^i^^f^yi S j 



4il \a+2u+l ^l — Q^Tn+l~2u 

A{m, I, a, u) 



(l + fi)2; Vl + 



.4(™,^a.u) = W(" + ':'V° + 'U"V'.'| ("4) 



where 



a + uj \a + u/ \u J \u 
and a = (ao, ai, . . . , an^i) € w'it/i a^- = — mj = kj — Ij. 

The proof of Proposition 6.10 is found in Appendix F. The respective special cases, like the 
four-dimensional Grosse-Wulkenhaar model, can easily be read off from this general expression. 



6.3 Power Counting 

The power counting theorem for general non-local matrix models was proven by Grosse and Wulken- 
haar in [GW05a]. In the matrix basis, every Feynman diagram of the duality covariant field theory 
is represented by a ribbon graph, whose topology is decisive for the question of whether or not it is 
divergent. Power counting in a dynamical matrix model depends crucially on this topological data. 
For a regular matrix model, the power counting degree of divergence for an iV-leg ribbon graph G 
of genus g with V vertices and B loops carrying external legs is given by [GW05a] 

Lo{G) = D + V {D - A) - I N {D - 2) - D {2g + B - 1) . (6.15) 

Here we briefly recall the role played by the asymptotic behaviour of the propagator in the derivation 
of this power counting theorem. For this, one uses multiscale analysis of the Schwinger parametric 
representation of the propagator, which works in both position and matrix space. 

The slicing of the propagator is defined as 

A = V A' through / ds = V / ds (6.16) 

with an arbitrary constant M > 1. This leads to a scale decomposition of the amplitude Ac of any 
given Feynman graph G as 

AG = J2Ah, (6.17) 

where / = {ii} runs through all assignments of positive integers ii to each line i of G. One then 
seeks appropriate bounds on the sliced propagators. 

For the i-th slice, the main bounds in matrix space are given by [RVTW06, GMRVT06] 

|^mr^;fc^| < M-^* g H^^+^+'^+'l'^ (6.18) 

5^ max I At,,^,, I < i^'M-^^e-'^'^^-'Mi (6.19) 

■'— ' n,k 
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for some positive constants K,K' and c, c', where we have introduced the £^-norm ||m||i := mg + 
mi + • • • + rrin-i- Perturbative power counting amounts to finding which summations cost a factor 
M^*^* through (6.18), 

-1 ]\/f2ni 

and which cost 0(1) due to the bound (6.19). Integrating out loops at higher scales of a graph 
then gives effective coupling constants in powers of M. The important point is that the faster 
the propagator decays, the smaller is the contribution of the integration over internal lines to 
effective coupling constants. This in turn reduces the number of divergent graphs which require 
renormalization. 

As an immediate application of this result in the present context, we can straightforwardly 
establish the power counting theorem for the 1+1-dimensional self-dual Grosse-Wulkenhaar model. 
In this case the matrix model is local with propagator A^^.^^ = 5mk ^ni A^'^\m, n) given by 

A^'^Vm,™) = — ^- : ^ (6.21) 



/•OO p\ 

I (|g g-2£;(m+n+l)s-ie-'«M^s = X ds (.-"^^ (m+n+l) s- i c ~ ■ 

Jo Jo 

,-1 



where -fC = i(e 2£;(m+n+i) le "-/x _ ^ Slicing this propagator as in (6.16), we easily find 
that the i-th slice can be bounded as 

|A('^)^(m,n)| < \K\M-^'{M^-1) ^-2E (m+n+l) M^^' _ (g_22) 

The Minkowski space propagator thus has the requisite exponential decays (6.18)-(6.19), and hence 
the perturbative multiscale renormalization in this case can be treated exactly as in the Euclidean 
setting [RVTW06]; we expect renormalizability to hold in this case. The case < 1 is much more 
difficult; in the Euclidean case the coupling il. fiows very rapidly to the self-dual point il. = 1, and 
it would be interesting to see if this is also the case for the hyperbolic self-dual point. 



6.4 Asymptotics 

As discussed in §6.3, the asymptotics of the propagators play an important role in perturbative 
renormalization. In this paper we are also interested in determining to what extent the complex 
matrix basis is applicable to the perturbative analysis of the duality covariant field theories; here 
the asymptotics also give crucial information. However, the asymptotic behaviours of the hyper- 
bolic parts of the propagators are difficult to investigate due to the oscillatory behaviours of the 
integrands. 

For example, consider the Grosse-Wulkenhaar model in four-dimensional Minkowski space with 
propagator given by Corollary 6.8. Introducing short variables = Xk — Uk and long variables 

= + Uk for A; = 0, 1, and using elementary hyperbolic and trigonometric identities, we can 
write this propagator in the form 

^'"'^-■^^ - -^^["^-''''^)^) («^^^' 

X exp J cot(s) \\uo\\u - J tan(s) ||i7o||m 
X exp ~^ coth(s) ||wi||e — tanh(s) ||i'i||e 
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where we set E = B for simplicity. The integral is sliced in the usual way to get 

^ ' ' (27r)2 sin(2s) sinh(2s) ^ ' 

X exp J cot(s) IIuoIIm - J tan(s) ||i7o||m 
X exp ~ ^ coth(s) ||mi||e — tanh(s) Hi'iHe 
with M > 1. 

The Euclidean part of the modulus of the integral (6.24) can be easily bounded from above 
by maximizing each of the hyperbolic functions in the integrand on the interval [M~'^^ , M~'^^^~^'^]. 
The factor e"? H^ilH takes its maximum at s = M'"^^ where tanh(s) = M'^-'' - ^ M'^^ + 

0(^{M~^^)^) < c' M~'^^ for some constant c' > 0, while e'^^- coth(s) ||iti||| ^g^j^gg maximum value 
at s = with coth(s) < M'^^^'^^ + M'^O-i) < c" M^J and some constant c" > 0. The 

function sinh(2s)~'^ can be bounded from above by M^-', and in this way one arrives at the very 
rough bound 

\A^''^^iu,v)\ < KM2^'e-'=(*^'^"-^llE+^^"'^ll"^llE) (6.25) 



Jm-2j sm(2s) V 4 4 



for some positive constants K and c. This reproduces the first required bound which gives expo- 
nential decay in both short and long variables in the Euclidean plane [GMRVT06]; in particular, 
integrating over long Euclidean coordinates costs a factor M^-' while short Euclidean coordinates 
cost M~'^K However, the asymptotic behaviour of the full propagator remains unclear; the hyper- 
bolic part of the integrand is oscillatory, so that more sophisticated methods are needed to bound 
this integral. 

There is a special case in which one can deduce the qualitative behaviour. The propagator of 
the critical, regularized massless LSZ model in 1-1-1 dimensions can be written using Proposition 6.1 
as 



where the integration contour has been rotated as s — )• s e Substituting 



u 



i-Ell^c - y|||(coth(2s^) - l) (6.27) 



we get 



A^^'Aix^y) = -^Q~-^--^y / dn 



' e---^«i^of|||^-y|li) , (6.28) 



47r V 2 

with Kq{^z) the modified Bessel function of the second kind of order 0. 

This implies that there is still a logarithmic ultraviolet divergence at x = y due to the singular 
behaviour of Kq{z) at z = 0. Since [AS70, 9.7.2] 



^o(^) = y^e-^ {l + 0{z~')) for 



oo , (6.29) 
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this also implies that the propagator A|^2'_j)(33, y) has an asymptotic exponential decay in the short 
variable u = x — y only for 

^R^ {\\u\\l) > , (6.30) 

and thus only for \'d\ < ^. We believe that for a < 1 the asymptotic exponential decay in the 
long variable v = x + y also persists as long as < ^. We conclude that the propagator has a 
worse behaviour in Minkowski space than in Euclidean space, but we can control its asymptotic 
behaviour with the help of the parameter i9. Regarding the restriction |??| < ^ as being part of 
the regularization of the field theory, one could then try to carry out the perturbative multiscale 
renormalization of the Minkowski space duality covariant field theory. 

In the matrix space representation there is a similar problem, since the integrand in (6.11) is 
oscillatory. Thus bounding the magnitude of the integral by an integral over the magnitude of 
the integrand possibly produces a big error and might lead to poor estimates of the asymptotic 
behaviour. One can use this approximation to show that the Minkowski space Grosse-Wulkenhaar 
propagator at = ^ has an exponential decay in each index separately, as in (6.18). To find 
the other bounds, however, one has to take care of the oscillating behaviour of the integrand. The 
asymptotics of the special case (6.28) for < ^ raises the hope that the propagators at hand may 
have such an asymptotic behaviour in position space so that the matrix basis is applicable.^ 
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A Proof of Proposition 3.18 



The complex Landau wavefunctions are built on tensor products of the complex harmonic oscillator 
wavefunctions fm"^^ as 



(A.l) 



/™"^(^) = V^w[|/(f''))(/f-)|] (x). 
Using (2.17) and (3.4) we get 

^iEkx/2 ^-lE^{it+k/2)^+(t~k/2f) 



E E, 



1 



47r V 27r V2"^+" ml n\ 
X H^{^/E^{t + k/2)) Hn{^/E^{t - k/2)) 



The generating function for the Hermite polynomials 

^-^ ml 



m=0 



(A.2) 



(A.3) 



^As these propagators are duality covariant, they also have a similar decay in momentum space. 
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is used to obtain the generating function for the complex matrix basis functions as 



V J—/ V 1 1 L t I L t 



m,n=0 

27r J 

X e "5 {t+fc/2)W-2r, (t-fc/2)) 



^ 2e5^'»(-<^- +2C-L"'+2^-V'-2^0 ^ (A.4) 

where we used the complex light-cone coordinates (3.20). The complex matrix basis functions can 
now be obtained by taking suitable derivatives with respect to the variables ^ and r/ to get 



m-\-n 



C=r,=0 



(m — p)! (n — p)!p! 



where we assumed m > n. This last sum can be identified with an associated Laguerre polynomial 

L^(.) = V . „ , i-zy (A.6) 



9=0 ('^-9)K^ + 9)! 



by shifting p ^ q = n — p. We finally arrive at 



TT V m! 



An identical calculation for n > m leads to the same result with + •(-)• — and m n, yielding (3.19). 



B One-Loop Effective Action in Matrix Space 

We will now reconstruct a classic result in quantum electrodynamics using the ??-regularization and 
the complex matrix basis. In his seminal paper [SchSl] Schwinger calculated the effective action 
for both a Dirac field and a Klein-Gordon field of charge e in a uniform external electromagnetic 
background in four spacetime dimensions. In a pure electric field E the one-loop correction for the 
Klein-Gordon theory (before charge renormalization) is given by the Lagrangian 

£W = ^ rd..-=^e-^^^f^^-l). (B.l) 
167r2 Jo Vsin(e£;s) J ^ ' 

By deforming the contour of integration above the real axis one picks up the poles at s = s„ = 
wk/cE for n G N by the residue theorem; this leads to the famous formula for the probability 
per unit time and unit volume 2 0'm(£^^^) to create a particle-antiparticle pair in the scalar field 
theory. We will now show that the regularized matrix basis approach leads to the same result quite 
effortlessly. We work throughout in the notations of §4. 



34 



The generating functional for connected graphs W[J, J*] is defined via the vacuum-to- vacuum 
amphtude (4.1) in the presence of the external sources J and J* as 

W[J,J*] = -i log Zo[J, r] . (B.2) 

It can be expressed in terms of the causal propagator (4.2) as [Sch51] 

W[J,r] = J d^x J d^y r{x)A,{x,y)J{y)- i\ogdet{Ap^A,) , (B.3) 

with Ap = Ac|^_Q the usual Feynman propagator. By using the ??-regularization we can write 

W[J,J*] = j d'^x j d'^y J* {x) M-'^^x, y) Jiy) - i log det ( ) ^ , (B.4) 



which is understood in the limit k — t- 0"*^ with 



(B.5) 



(B.7) 



K^-H^ Jk ei'^PSd -k) - e-i^/x2 + (iax)2 
The effective action is now defined as the Legendre transformation 

r[,^ei,</'d] = W[J,J*] - J d^x J{x)4>l,{x) - j d^x r{x)(t>,i{x) (B.6) 
of W[J, J*] with respect to the "classical" fields (l)c\{x) and (l)*^i{x) defined by 

<^ci(^) = ^-^^ = j A('=.^)(^,y) J(y) , 

^:i(-) = ^^g^ = / d^y J*(y)AM)(y,.). 

These equations may be inverted to give 

J{x) = {Kl-^?)^ct>,,{x) and J* {x) = -{Kf^ - i?),<Pl,{x) , (B. 
and inserting this into (B.6) yields 

r[</.,l,0:i] = - / j d^y{\<'^-^?)^<l^:,{x)A^-^^\x,y){Kl-^i^)My) 
i log det ^ 



d'^x {{Kl-^?)Mx))4>Ux) + j d^x {{Kl-i?)^4>li{x))4>ci{x) 
?o['/'ci,</'d]+ ilogdetf ^^-^1^ ) . (B.9) 



This is the juU effective action of the quantum field theory; the quantum mechanical content is 
completely captured by the one-loop correction 

ilogdetf J =T^[0,0]. (B.IO) 



-5,2 — fl"^ / K 
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We define W[0,0] =: J d^a; C^^\x), with the one-loop effective Lagrangian. The probabihty 
that no pair gets produced out of the vacuum is given by |(0,out|0,in)| = 

The effective action is given by 



-2am{iy[o,o]) 



W[0,0]= iTrlog(-|^^j^ , (B.ll) 

and the eigenvalue equation for the operator (K^— /i^)^ is given by (4.17). We adhere to Schwinger's 
convention by substituting E — t- eE/2. Using the identity 



which is valid for Jm{a) > and Jm{b) > 0, the effective Lagrangian can be obtained through 
C^'\x) = i(a^|log( \x) (B.13) 



oo Ac, r a2 



m,n=0 ^ ' 

The integration over parallel momenta gives 

We can now use Lemma 4.13 with x = y and a = 1 to obtain 



1 ds _:„„2 / eE 1 



IS 



IGvr^ Jq Vsinh(e-Es) s 



(B.14) 



which is independent of x. The integral converges near infinity since fif. has a small imaginary 
part, and near due to the ^ subtraction of the free scalar propagator. By rotating the integration 
contour as s — t- — is, and taking the limit k — t- 0"*", this Lagrangian coincides with Schwinger's result 
(B.l). The case of four-dimensional Dirac fields can be treated in the same way, by transforming 
the spinor propagator to the scalar propagator; see [Fisll, App. F] for details. This analysis again 
exemplifies the fact that the matrix basis provides an easy way of doing otherwise cumbersome 
calculations in quantum electrodynamics. 

C Proof of Lemma 4.13 

For m > n the explicit expression for the first eigenfunctions on the left-hand side of (4.14) is 

/f/)(-) = (-^r^l^^^e-^^^^'^-'^' [^,x^^Y-LT-{E,.f .^^) , (C.l) 
while the second eigenfunctions have a similar representation 
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with the notations (2.31) and (3.20). These representations can also be used for n > m due to the 
identity 

(-1)" r'"-" L™-"(r2) = (-1)"^ r"-"^ ^ L'lr"'(r'^) . (C.3) 

n! 

The sum over n thus has the form 

E/S'w /£«'(«)»" = ^ (^■"-'-'5')'% -.,,.r..'>.,r»'-")/^ (c,4) 

vr m! 

n=0 



X 

n=0 



oo 



It can be done explicitly by using the identity [Han75, 48.23.11] 

oo 

n\ e Lr "(6 l^'iri^l) = fc! e ^^'^ (1 - c)™"^^ ' ((1 " ^ c) (77 c - l)/c) (C.5) 



n=0 



With = m, e = i?. x^^\ r? = i^^ and c = a/E, x^^^ y V \ This yields 

00 p 

E ^(^) ^(2/) = - e e'=«'' a- (7? + e - c^,? - c-i) , (C.6) 

n.=0 

which after some elementary algebra gives (4.14). 

D Matrix Basis in Momentum Space 

The complex Landau wavefunctions have special symmetries which will be useful in analysing their 
Fourier transforms. 

Proposition D.l. The complex Landau wavefunctions satisfy the relations 

f!^^\E-H,E-'x) = Ef^l''-'\t,x) , (D.2) 
fL't\-t,x) = (-l)™-"/(i-)(t,x) , (D.3) 
/(f„'')(i,-x) = fig^\t,x) , (D.4) 

f!nn'H^,t) = (- i i^-"^ (t, x) . (D.5) 

Proof. The relation (D.2) follows directly from the explicit expression (3.19) by noting that E and 
x^^^ occur only in the combinations \/^ and E xf^^ . Time reversal t ^ —t only affects the 
term involving x^^^ , ^ — t- —x"^\ , = —x^^^ , which gives (D.3). Parity x — )• —x sends 

° — sgn(m— n) sgn(m— n) — sgn(n— m) ' >=> \ j J 

x^^^ , ^ — )• x^^^ I N = x^^^ , ^, which shows (D.4). Under interchange of t and x, we find 

— sgn(m— n) sgn(m— n) — sgn(n— m)' ^ ' ° ' 

4^) ^±ie-i'^x^-''\andthusV^xf ^ v^l^i^"''^) and i?^ xf xL^^ ^ 4"'') xL"''^ 
Putting these transformations into (3.19) proves (D.5). □ 

Proposition D.6. The Fourier transformation of the complex Landau wavefunction fmn (x) is 
given by 

( • \m—n 

nf'^r!^r?^m = fi]/^^\k) = ^-^^ fLl^\E''k) (D.7) 

with E-^k = (A;\A;°). 
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Proof. Denote momentum space derivatives as := Using the explicit forms of the hyperbohc 
and Euchdean space wave operators given in §2.1, in Fourier space we find that these operators 
have the form 

i- I d2^(p2<^)(^)e-'^- = h I (D-8) 

= ((fco -kl)+2iE {k^ 3^ - fci 5°) + {di - dl))j^[cP]{k) 

and 

^ J d2a:(Pf<^)(a;)e-'^- = ^ J (D.9) 

= {iki + kl) + 2[Eik''d^ + k^d'^)-E''{di + df))Tmk) . 

From the expHcit forms of the regularized wave operators (2.26), this gives 

^[P2(??)0](fe) = e'^E^-{e'^di+ e''^dj) + 2iE-^{e'^k^d^+ e-'^k^d^) 

+ ie-'''ki+ e''^kl))T[(l)]ik) 

= e'^'^ E^V^i-^mjik) , (D.IO) 
where the differential operator 'P^(— i?) has the same form as P^(— i?) with the substitutions — ?• d^, 

( F ) 

— )• /c^ and E — t- E . On the other hand, by substituting (j) = fmn we find 

:F[P\^)f'^r?A (fc) = ^E, [m + 1) ^[/(f^")] (fc) (D.ll) 

and thus 

^2(_^)^[/(f„.)] (fc) = AE-' (m + i) ^[/(f,-)] (fc) . (D.12) 

By Parseval's theorem the Fourier transforms of the matrix basis functions have the same normal- 
ization as the position space wavefunctions, from which we conclude the first equality of (D.7). The 
second equality of (D.7) follows from the symmetry relations (D.2)-(D.5). □ 

As a simple application of this result, we can establish that the Feynman propagator for the 
free Klein-Gordon theory in the complex matrix basis possesses the same mass-shell singularities 
as in momentum space. 

Lemma D.13. The Feynman propagator in the complex matrix basis is given by 
with 

and e ^ /or t9 = I - K > 0. 

Proof. We simply relate the Klein-Gordon operator in the different basis sets. One has 

n,m,k,l 
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and thus 

Gl:U = /d'- ld'yf^rrH-){-\{dl + f^')Jy)flf''\y) 

= I ^ 1^'^ I d'yer:\^)e-^'-{-\\k\\l + ,.l)e-^''-yfif^\y) 

= J d'k ^[/(^-^)] (fc)* ( - llfcll^ + fil) -F[/(f-)] {k) . (D.17) 

It follows that the Fourier transforms of the functions f!wn^ diagonalize C^^.^^, and the result now 
follows from (D.7). □ 

E Proof of Proposition 6.1 

The propagator is given by 
A('^''^)(a;,y) = {x\{ae'''K'^{^) + de'''K^{'&)- e-'^'fi^y^ly) (E.l) 



where i} = ^ — k > and we have set a = 1 — a. The (regularized) wave operators have the 
eigenvalue equations 

(aP2(z?) + aP2(^9))/i^-„)(a;o)=4ii;^(amo + ano + i)/itn^^ , 

(E.2) 

{a (P2)fc+i + a (Pf)fc+i)/iti(^fc) = 4i?fc {amk+ank + i) fr[f±{xk) , 

with fm^nki^k) the usual Landau wavefunctions and Bj^ G ]R_|_ for k = 1, ... ,n — 1. Using the 
identity 



f oo 

-1 _ / A „ „—sa 



ds e-'" (E.3) 







which is valid for 9^e(o) > 0, we find 



oo 



A('^'-)(a;,y) = "ie-'W e"^^'^ /S(^o) /^(^/o) e"^^^-^ (-'"»+^""+^) 

7Tio,no=0 



n— 1 oo 



xn( E /itnl(^^0/S^l(y;^Oe-^^^'=('^™'=+^"'^+^)) . (E.4) 

fc=l mfe,nfc=0 

By Lemma 4.13 the sum over uq gives 

- I ||a=o - yolli + (e-^^^--"^ -1)E {xo,yo)^ - e-^^^-"^ ixo ■ Eyo) 
lI^^ {E\\xq - yoWl - Asinh^ils a) E {xo,yo)^ - 2 sinh(4s ct) ixq ■ Eyo) . 



mo=0 



X 
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The sum over mo can be performed by using the identity (4.20) with t = e~'^''^-« to get 

^^P - o ■ uVo J\ E\\x^-y4l (E.6) 



27r sinh(2s£;_^) V 2 sinh(2s 

+ ( - e----^+ e---^^^^i^^^^) i.o-^yo 
\ smh(2s£/_^j / 

By using elementary hyperbohc identities, the terms proportional to (a;o, Z/o)i? can be simplified to 

^-'isE_,,a o^-2si;_., sinh^(2g£^_^a) _ cosh(25£^_^) cosh(2g£^_^) 

sinh(2s£;_^) sinh(2s£;_tf) sinh(2s£;_^) ^ ' ' 

where we defined i?-.,? := (1 — 2a')£'„,9 = {2(t—1)E_^. Likewise, the terms proportional to ixQ-Eyo 
can be rearranged to 

^-4sE_„a . ^-2sj;_., sinh(4g£;_^ff) ^ sinh(2g^_^) 

sinh(2s£;_.^) sinh(2s£;_^) ' ^ ' ^ 

The triangle relation ||a;o — yo||.| = ll^^olll + l|yo||| — 2(a;o, Vo)^ allows us to combine further terms. 
The sums over and for k = 1, . . . ,?i — 1 are treated in exactly the same way, and putting 
everything together we finally get (6.2). 



F Proof of Proposition 6.10 



The 2n-dimensional regularized LSZ wave operator in the matrix basis is given by (3.60)-(3.61) 
with 6j = 9 and 'D^^^.j.i '■= T^mri^ki for j = 0, 1, . . . , n — 1. Each of these operators have non-vanishing 
matrix elements only for 

'^i ~ "^i ~ ~ h ~- ^3 fo^ J = 0, 1, . . . , n — 1 . (F.l) 

This is due to the 50(1, 1) x S0{2)'^~^ symmetry of the action. We can thus eliminate n components 
and write instead 

n-l 

i=l 

with Q G Z". 

The n components of the wave operator (F.2) are independent and its eigenvectors are therefore 
products of the eigenvectors of the individual matrices. The mass term is already diagonal and so 
are the terms proportional to 0,. Thus for every a G Z we seek solutions of the eigenvalue equations 

oo 

Sr^ j^{l/2) tM) _ rjia) (-p o) 

^m,m+a;l+a,l'-^lv — ^ '-^ mv ■ \^ 

1=0 

This equation has been solved in [GW05b]. The eigenvectors are given by 

An ' 



(„) I fa + m\ fa + y\ / 2VUy+^ fl-ny'-^y ^(-m,-y 
^'m^ = ^^^ m )\ y 1 + a 



(F.4) 
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and the eigenvalues are 



V 



(2y + a + l) (F.5) 



for y E No- As expected, this is the usual harmonic oscillator spectrum. The hyper geometric func- 
tion 2-Fi appearing in (F.4) with negative integer values in its first two arguments is an orthogonal 
Meixner polynomial. In particular, Umv is symmetric in its lower indices. 

For the full wave matrix the addition of the fi-term modifies the eigenvalues u — )• f ' with 

v' = ^(2y + 2aa + l) . (F.6) 
The complete matrix space wave operator va. D = 2n dimensions has the representation 



m,m+a 



n-1 



i=l 

where 



ra-l 

Umv = JJ UrrijVj (F. 



j=0 

and 

n—l n—1 
i=l 1=1 

with Uj G Nq. From the orthogonality relations for the Meixner polynomials it follows that 

In the following we will use the notation Umv = Um\y) where v and y are related by (F.5). 
Using the Schwinger parametrization this yields the propagator 

oo n—l -. n—l 

= _ei^— /"^ dt e^*''"'('^"°+5)-*Er=/(^°«+|)-*^ fFll) 

m Jo ^ ■ ^ 

oo n—l oo 

x( ^ e'*^'^^o^(-)(yo)<"°^(2/o)) n (E 

yo=o i=i j/i=o 

The sum over yo can be performed by using the explicit formula for the eigenvectors (F.4), and the 
hypergeometric identity [GW05b] 



oo 

^ I rv -\- H \ I — m — II 

W I 2F1 



Oi + y\ „ ( -m, -y 



1 2^1 I 1 I 
y J \ 1 + a 



y=Q 

^ {l-{l-w)zr^' ( -m,-l 

(1 _ ^)Q+m+«+l 2 1 I 



-h-y 

l + a 



(F.12) 



9 

z w 



{l-{l-w)z)\ 



for Izl < 1 
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with z 



it e ' 



(1 - nf (1 + ny^ and w = -An (1 - ny 



After some algebra this leads to 

oo 

yo=o 



(F.13) 



X2F1 



mo 



-mo, -lo 
1 + ao 
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In 



2 / l + 0\2 



1 - n 



it c 



1-e 



it e ■ 



Now we substitute s = e * (with Jacobian s ^) and use the expansion of the hypergeometric 
functions 



2F1 



-m, —I 
1 + a 



min(m,i) 

E 

tt=max(0,— a) 



ml U a\ 



(m — u)\ {I — u)\ {a + u)\ u\ 



(F.14) 



After a bit of algebra the various factorial terms can be recombined into the quantity (6.14), and 
we find 



(F.15) 



2/0=0 



min(mo,io) 



E 



i e "=-uo ^-|^ _ g-i e 



ito=max(0,— do) I 1 



(i+n) 



ao+"iO+'o+l 
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(1 + 0)2 



ao+2no+l / 1 — O \ mo+lQ-2uo 



1 + 



A{mo,lo,ao,uo) 



The sums over yi for i = 1, . . . , n — 1 are performed in a completely analogous way. The only 
difference between the Euclidean and hyperbolic parts of the propagator is the additional factor 
— i e ''^ in the exponential of yo- This simply changes it e — t- — t and ^ — )• s everywhere in 
the above derivation, and we arrive finally at the expression (6.11). 



References 

[AS70] M. Abramowitz and I. Stegun. Handbook of Mathematical Functions. Dover Publica- 

tions, New York, 1970. 

[Bah04] D. Bahns. Perturbative methods on the noncommutative Minkowski space. PhD 
Thesis, 2004. DESY-THESIS-2004-004. 

[BahlOa] D. Bahns. Schwinger functions in noncommutative quantum field theory. Ann. Henri 
Poincare, 11:1273-1283, 2010. 

[BahlOb] D. Bahns. The ultraviolet infrared mixing problem on the noncommutative Moyal 
space. Preprint, 2010. arXiv: 1012.3707 [hep-th]. 

[BDFP02] D. Bahns, S. Doplicher, K. Fredenhagen, and G. Piacitelli. On the unitarity problem 
in space/time noncommutative theories. Phys. Lett. B, 533:178-181, 2002. 



42 



[BFS85] LA. Batalin, E.S. Fradkin, and Sh.M. Shvartsman. Quantum electrodynamics in 
external constant field. Nucl. Phys. B, 258:435-467, 1985. 



[BGM+07] A. P. Balachandran, T.R. Govindarajan, G. Mangano, A. Pinzul, B.A. Qureshi, and 
S. Vaidya. Statistics and UV/IR mixing with twisted Poincare invariance. Phys. Rev. 
D, 75:045009, 2007. 

[CR77] P. Candelas and D.J. Raine. Feynman propagator in curved spacetime. Phys. Rev. 

D, 15:1494-1500, 1977. 

[Dav99] E.B. Davies. Pseudo-spectra, the harmonic oscillator and complex resonances. Proc. 
Roy. Soc. London Ser. A, 455:585-599, 1999. 

[DeW75] B.S. DeWitt. Quantum field theory in curved space. AIP Conf. Proc, 23:660-688, 
1975. 

[DFR95] S. Doplicher, K. Fredenhagen, and J.E. Roberts. The quantum structure of spacetime 
at the Planck scale and quantum fields. Commun. Math. Phys., 172:187-220, 1995. 

[DK04] E.B. Davies and A.B.J. Kuijlaars. Spectral asymptotics of the non-self-adjoint har- 

monic oscillator. J. London Math. Soc, 70:420-426, 2004. 

[Dun09] G.V. Dunne. New strong field QED effects at ELI: Nonperturbative vacuum pair 
production. Eur. Phys. J. D, 55:327-340, 2009. 

[FG81] E.S. Fradkin and D.M. Gitman. Furry picture for quantum electrodynamics with pair 

creating external field. Fortsch. Phys., 29:381-412, 1981. 

[FGS91] E.S. Fradkin, D.M. Gitman, and Sh.M. Shvartsman. Quantum Electrodynamics with 
Unstable Vacuum. Springer- Verlag, 1991. 

[Fisll] A. Fischer. Matrix model formulation of LS-covariant noncommutative quantum field 

theories on Minkowski spacetime. PhD Thesis, 2011. Leibniz Universitat Hannover. 

[FS09] A. Fischer and R.J. Szabo. Duality covariant quantum field theory on noncommutative 

Minkowski space. J. High Energy Phys., 02:031, 2009. 

[FSll] A. Fischer and R.J. Szabo. UV/IR duality in noncommutative quantum field theory. 

Gen. Relativ. Gravit., 43:2509-2522, 2011. 

[Git77] D.M. Gitman. Processes of arbitrary order in quantum electrodynamics with a pair 

creating external field. J. Phys. A, 10:2007-2020, 1977. 

[GL07] H. Grosse and G. Lechner. Wedge-local quantum fields and noncommutative 

Minkowski space. J. High Energy Phys., 11:012, 2007. 

[GL08] H. Grosse and G. Lechner. Noncommutative deformations of Wightman quantum 

field theories. J. High Energy Phys., 09:131, 2008. 

[GMOO] J. Gomis and T. Mehen. Space/time noncommutative field theories and unitarity. 
Nucl. Phys. B, 591:265-276, 2000. 

[GMRVT06] R. Gurau, J. Magnen, V. Rivasseau, and F. Vignes-Tourneret. Renormalization of 
noncommutative (/>| field theory in x space. Gommun. Math. Phys., 267:515-542, 
2006. 



43 



[GRVT06] R. Gurau, V. Rivasseau, and F. Vignes-Tourneret. Propagators for noncommutative 
field tlieories. Ann. Henri Poincare, 7: f 601-1628, 2006. 

[GW03] H. Grosse and R. Wulkenhaar. Renormalization of theory on noncommutative 
in the matrix base. J. High Energy Phys., 12:019, 2003. 

[GW05a] H. Grosse and R. Wulkenhaar. Power-counting theorem for non-local matrix models 
and renormalisation. Commun. Math. Phys., 254:91-127, 2005. 

[GW05b] H. Grosse and R. Wulkenhaar. Renormalization of (f)^ theory on noncommutative 
in the matrix base. Commun. Math. Phys., 256:305-374, 2005. 

[Han75] E. Hanson. A Table of Series and Products. Prentice-Hall, Englewood Cliffs, N.J., 
1975. 

[HI09] T. Heinzl and A. Ilderton. Exploring high-intensity QED at ELI. Eur. Phys. J. D, 

55:359-364, 2009. 

[ILMIO] A. Ilderton, J. Lundin, and M. Marklund. Strong field, noncommutative QED. 
SIGMA, 6:041, 2010. 

[LS02] E. Langmann and R.J. Szabo. Duality in scalar field theory on noncommutative phase 

spaces. Phys. Lett. B, 533:168-177, 2002. 

[LSZ03] E. Langmann, R.J. Szabo, and K. Zarembo. Exact solution of noncommutative field 
theory in background magnetic fields. Phys. Lett. B, 569:95-101, 2003. 

[LSZ04] E. Langmann, R.J. Szabo, and K. Zarembo. Exact solution of quantum field theory 
on noncommutative phase spaces. J. High Energy Phys., 01:017, 2004. 

[MVRSOO] S. Minwalla, M. Van Raamsdonk, and N. Seiberg. Noncommutative perturbative 
dynamics. J. High Energy Phys., 02:020, 2000. 

[PialO] G. Piacitelli. Quantum spacetime: A disambiguation. SIGMA, 6:073, 2010. 

[Rin03] A. Ringwald. Fundamental physics at an X-ray free electron laser. Erice 2001: Elec- 
tromagnetic Probes of Fundamental Physics, pages 63-74, 2003. Eds. W. Marciano 
and S. White, World Scientific, River Edge. 

[Rit70] V.I. Ritus. Mass separator and exact Green's function of the electron in an intense 

field. Pisma Zh. Eksp. Teor. Fiz., 12:416-418, 1970. 

[Rit78] V.I. Ritus. Eigenfunction method and mass operator in quantum electrodynamics of 

a constant field. Sov. Phys. JETP, 48:788-799, 1978. 

[Riv07] V. Rivasseau. Constructive matrix theory. J. High Energy Phys., 09:008, 2007. 

[RVTW06] V. Rivasseau, F. Vignes-Tourneret, and R. Wulkenhaar. Renormalization of non- 
commutative (;/)^-theory by multi-scale analysis. Gommun. Math. Phys., 262:565-594, 
2006. 

[Sch51] J. Schwinger. On gauge invariance and vacuum polarization. Phys. Rev., 82:664-679, 

1951. 

[Sol07a] M.A. Soloviev. Noncommutativity and 6'-locality. J. Phys. A, 40:14593-14604, 2007. 



44 



[Sol07b] M.A. Soloviev. Star product algebras of test functions. Theor. Math. Phys., 153:1351- 
1363, 2007. 

[Solos] M.A. Soloviev. On the failure of microcausality in noncommutative field theories. 

Phys. Rev. D, 77:125013, 2008. 

[Sza03] R.J. Szabo. Quantum field theory on noncommutative spaces. Phys. Rept., 378:207- 

299, 2003. 

[Wanll] Z. Wang. Construction of 2-dimensional Grosse-Wulkenhaar model. Preprint, 2011. 
arXiv:1104.3750 [math-ph]. 

[Zahlla] J. Zahn. Divergences in quantum field theory on the noncommutative two-dimensional 
Minkowski space with Grosse-Wulkenhaar potential. Ann. Henri Poincare, 12:777- 
804, 2011. 

[Zahllb] J. Zahn. Ultraviolet-infrared mixing on the noncommutative Minkowski space in the 
Yang-Feldman formalism. Preprint, 2011. arXiv: 1105. 6240 [hep-th]. 



45 



